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ABSTRACT 


In  1967  and  1968  instrumented  towers  and  thermograph  networks 
were  placed  in  operation  in  the  cities  of  Edmonton  and  Calgary,  Alberta, 
for  continuous  measurement  of  winds  and  temperatures  in  prairie  city 
environments.  In  addition,  a  mobile  unit  was  provided  for  horizontal 
temperature  traverses  and  tape  samplers  yielded  continuous  horizontal 
of  smoke  concentrations  at  several  points  in  each  city.  Some  of  the 
data  from  these  sources  are  used  in  analyses  of  the  structure  of  the 
urban  heat  island  in  Edmonton  and  Calgary  and  in  analyses  of  the  con¬ 
tribution  of  selected  meteorological  parameters  to  the  total  variability 
of  daily  smoke  concentrations  in  Edmonton. 

A  statistical  model  is  derived  for  the  prediction  of  mean  daily 
output  of  smoke  in  Edmonton.  However,  the  derived  discriminant  functions 
fail  and  this  failure  is  attributed  to  either  (a)  improper  choice  of 
predictors  or,  (b)  poor  estimates  of  the  predictand.  Seasonal  temper¬ 
ature  variations  do  not  account  for  a  significant  portion  of  the 
variation  in  smoke  output.  A  possible  explanation  for  this  is  that 
space  heating,  which  is  temperature  dependent,  is  not  a  major  source 
of  smoke,  due  to  the  use  of  natural  gas. 

Using  wind  speed,  horizontal  temperature  differences,  and 
vertical  temperature  differences,  a  linear  multiple  regression  equa¬ 
tion  is  derived  for  the  prediction  of  mean  daily  smoke  concentrations 
in  the  city  of  Edmonton.  The  equation  explains  about  407>  of  the  total 
variation  in  smoke  concentration,  and  shows  considerable  stability  in 
a  test  on  independent  data. 


Sip;. 


A  significant  change  in  temperature  is  found  as  one  approaches 
the  center  of  the  city  from  the  outskirts.  This  difference  is  attribu¬ 
ted  to  differences  in  heat-storage  capacity,  differences  between  in¬ 
coming  and  outgoing  radiation  due  to  pollution,  and  the  excess  heat 
due  to  human  activities  in  the  city.  The  urban  heat  island  is  divided 
into  at  least  two  parts  by  the  valley  of  the  North  Saskatchewan  River. 

Analysis  of  the  vertical  temperature  distribution  in  Edmonton 
reveals  that  near-adiabatic  conditions  prevail  during  the  daylight 
hours  in  the  summer;  and  isothermal,  or  a  weak  positive  or  negative 
lapse  prevails  throughout  the  rest  of  the  year.  Temperature  inversions 
of  moderate  strengths  occur  throughout  the  year  at  night  and  in  the 
early  morning  hours. 

Analysis  of  the  vertical  temperature  structure  in  Calgary 
shows  that  adiabatic  to  superadiabatic  conditions  prevailed  during  the 
daylight  hours  in  March,  1968,  through  September,  1968.  Moderate  to 
strong  inversions  occur  throughout  these  months  at  night  and  in  the  early 
morning  hours. 

When  results  from  other  sites  are  studied,  it  is  evident  that 
the  Edmonton  Tower,  more  so  than  the  Calgary  Tower,  represents  a  typical 
city  environment  possibly  because  of  the  peripheral  location  of  the 


latter . 
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INTRODUCTION 


Pollution  History 

The  quality  of  the  atmosphere  on  which  existing  terrestial 
forms  of  life  are  dependent  has  been  recognized  as  important  only 
during  the  past  few  decades.  It  can  be  supposed,  however,  that  smoke 
from  volcanoes,  forest  fires,  and  domestic  heating  and  cooking  arrange¬ 
ments  were  locally  troublesome  to  our  human  ancestors.  But  it  is  very 
unlikely  that  such  circumstances  were  more  than  incidents,  until  family 
units  began  to  live  in  communities;  only  then  could  man  produce  waste 
enough  to  affect  the  neighborhood.  Not  much  is  known  from  the  Medieval 
period,  but  from  accumulated  knowledge  of  domestic  heating  practices 
and  limited  industrial  activities  one  can  draw  the  conclusion  that  air 
pollution  was  not  a  serious  problem  in  the  villages  and  towns  of  the 
time.  The  problem  of  pollution  was  first  recognized  following  the 
replacement  of  wood  by  coal  as  the  prime,  source  of  energy.  Various 
limitations  and  prohibitions  to  the  use  of  coal  were  proclaimed  official¬ 
ly  and  there  are  even  records  on  punishment  of  offenders.  Until  the 
early  part  of  the  twentieth  century  coal  smoke  was  the  major  contributor 
to  air  pollution  and  even  now  coal  is  still  used  extensively  in  many 
industrialized  areas  of  the  world.  The  use  of  petroleum  products  in 
the  forms  of  gasoline  and  oil  has  been  enormously  accelerated,  due  to 
the  internal  combustion  engine.  The  combustion  of  oil  and  gas  has 
diminished  the  coal  smoke  nuisance  but  has  instead  produced  "new"  kinds 
of  air  pollution.  One  well-known  example  is  the  "Los  Angeles  smog". 

This  results  from  the  exposure  to  sunlight  of  products  of  petroleum 


mixed  with  oxides  of  nitrogen.  Although  most  of  the  pollutants  related 
to  petroleum  processing  and  use  have  toxic  and  irritative  properties  of 
low  order,  their  photochemical  reaction  products  even  at  very  low  concen¬ 
trations  may  affect  life.  A  new  phase  of  the  air  pollution  problem  will 
arise  with  the  introduction  of  energy  sources  such  as  nuclear  power  and 
solar  energy.  No  one  knows  yet  for  certain  how  and  to  what  extent 
the  radioactive  by-products  of  nuclear  power  will  affect  the  biological 
environment . 

The  effects  of  air  pollution  can  be  manifold.  It  is  well  known 
that  air  pollution  will  reduce  visibility  due  to  the  scattering  of  light 
from  the  surfaces  of  the  airborne  particles.  The  degree  of  reduction 
is  related  to  the  optical  airmass,  particle  size  and  the  density  of  the 
medium.  The  visibility  is  so  poor  at  times  in  many  industrialized 
cities  that  transport  and  other  municipal  activities  have  to  be  curtail¬ 
ed.  Material  damage  is  another  widespread  effect  of  air  pollution, 
the  value  of  which  has  been  estimated  to  billions  of  dollars  annually 
in  the  United  States  alone1.  A  number  of  plants  and  crops  are  very 
sensitive  to  certain  pollutants  and  if  good  air  quality  is  not  provided 
for,  this  can  result,  for  example,  in  decreased  size  and  yield  of 
fruits  and  destruction  of  flowers.  Under  extreme  conditions  the  air 
can  become  polluted  to  the  extent  that  it  is  deadly.  Typical  examples 
from  the  1940's  and  1950 's  are  London,  Donora  and  Meuse  Valley.  Long, 
continued  exposure  to  light  or  moderate  concentrations  of  pollutants 

^R.E .  Munn,  "The  Interpretation  of  Air  Pollution  Data,  with 
Examples  from  Vancouver",  Canada,  Dept,  of  Transport.,  Met.  Branch, 
Technical  Circular  Series,  Cir-3454,  Tec-351,  1961,  p.  1. 
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may  have  physiological  effects  although  direct  proof  is  difficult  to 
establish.  The  many  cases  of  chronic  bronchitis  in  British  cities  and 
irritated  eyes  in  Los  Angeles  seem  to  be  closely  associated  with  air 
pollution.  Although  awareness  of  the  problem  of' city  air  pollution  is 
increasing  it  is  viewed  largely  as  a -problem  of  the  future.  Only  a  few  large 
cities  are  now  using  the  supply  of  clean  air  faster  than  nature  can  restore 
it.  Such  overuse  is  likely  to  occur  with  increasing  frequency  as 
populations  increase  and  migration  to  urban  areas  continues. 

Diffusion  Theory 

Serious  studies  on  turbulence  and  associated  diffusion  processes 

started  with  the  simple  experiments  on  the  flow1  of  w?ater  in  long  straight 

2 

pipes  described  by  0.  Reynolds  as  early  as  1883  .  According  to  Reynolds, 
turbulent  flow  could  be  modelled  as  a  simple  mean  motion  on  which  was 
superimposed  a  complicated  eddy  motion  of  an  oscillatory  character. 

In  his  experiments  he  managed  to  show  that  the  motion  became  turbulent 
when  the  so-called  Reynolds  number  exceeded  certain  critical  values. 

Much  of  the  early  work  on  atmospheric  turbulence  was  based  on  theory 
that  was  reminiscent  of  the  theory  of  viscosity  developed  in  the  kinetic 
theory  of  gases.  Like  molecules  the  eddies  transfer  momentum  from  one 
layer  to  another  by  breaking  away  from  the  main  stream  in  a  region  of 
high  velocity,  conserving  some  or  all  of  its  momentum  until  it  mixes 
again  with  the  mean  motion  in  a  region  of  lower  velocity.  Such  momentum 
transfer  would  show  up  as  an  increase  in  velocity  at  the  new  level  and 

_ 

O.G.  Sutton,  Micrometeorology,  A  Study  of  Physical  Processes 
in  the  lowest  Layers  of  the  Earth's  Atmosphere,  McGraw-Hill,  New 
York,  1953,  p.  56. 
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would  have  the  effect  of  eliminating  large  velocity  differences. 

This  theory  was  further  developed  by  Wilhelm  Schmidt  who  introduced  the 

coefficient  of  viscosity  which  he  called  "Austausch  K.oef f izient" 

("exchange  coefficient")  .  The  analogy  between  molecular  and  eddy 

motion  was  further  extended  by  Ludwig  Prandtl  in  1925,  who  brought  in 

the  mixing -length  hypothesis  by  introducing  a  quantity  L  which  is  a 

unique  length  characterizing  the  local  intensity  of  turbulent  mixing 

4 

at  any  level  and  which  is  a  function  of  position  .  Though  much  fruit¬ 
ful  work  has  been  done  using  the  preceding  hypothesis  it  became  evident 
through  recent  experimental  work  that  the  mixing- length  theory  had 
serious  defects.  In  many  turbulent  exchange  problems  it  has  been  aban¬ 
doned  in  favor  of  the  so-called  statistical  theories.  This  new  theory 
was  first  suggested  by  G.I.  Taylor”*  (1921.)  in  a  paper  on  the  problem 
of  diffusion  written  in  a  journal  devoted  mainly  to  pure  mathematics. 
The  theory  depends  largely  upon  the  concept  of  correlation  in  statis¬ 
tics,  which  is  an  attempt  to  measure  closeness  of  relationship  numer¬ 
ically.  In  his  work  from  1921  Taylor  studied  the  extension  of  the 
random  walk  problem  in  which  there  is  a  variable  correlation  between 
the  motion  of  a  particle  at  one  time  and  at  some  instant  later.  His 
famous  correlation  equation  is  now  the  starting  point  in  many  fluid 
mechanics  problems.  Further  advances  in  the  field  of  turbulence  were 


3 

Ibid. ,  p .  68 . 

^Ibid . ,  p .  72. 

-*G.I.  Taylor,  Proceedings  of  the  London  Mathematical  Society, 
Ser.  2,  Vol .  XX,  1921,  p.  196-212. 
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made  with  the  introduction  of  Monin- Obukhov ' s  similarity  theory  in 

1954  .  This  theory  states  that  there  exist  near  the  ground  a  velocity 

U* ,  a  length  L  and  a  temperature  T*  that  are  essentially  invariant  with 

height.  When  temperature,  wind,  and  height  are  expressed  nondimensionally 

as  functions  of  these  quantities,  a  series  of  nondimensional  equations 

results  that  is  of  very  general  validity  in  the  surface  boundary  layer. 

One  of  the  problems  of  the  similarity  theory  is  to  relate  various 

dimensionless  variables  to  the  dimensionless  height  ratio  z/L  where  L 

is  the  above  mentioned  length  (scaling  length). 

As  the  theory  was  not  fully  tested  by  Monin  and  Obukhov  in 

the  field,  Taylor  (1960)  examined  the  similarity  theory  using  data 

from  work  done  by  Rider  (1954)  and  Swinbank  (1955)  and  concluded  that , 

except  in  strong  lapse  conditions  where  some  doubt  remained,  the  simi- 

7 

larity  theory  of  Monin  and  Obukhov  applied. 

A  remaining  difficulty  is  the  question  of  relative  size  of 
the  eddy  conductivity  and  the  eddy  viscosity  EL.  which  is  import¬ 

ant  in  the  determination  of  the  scaling  length  L.  Evidence  and  argu¬ 
ments  have  occasionally  been  put  forth  indicating  that  the  ratio  K^/Kqj 

8 

varies  with  stability  (Swinbank  1955)  .  However,  other  investigators 

^J.L.  Lumley  and  H.A.  Panofsky,  The  Structure  of  Atmospheric 
Turbulence ,  John  Wiley  &  Sons,  New  York,  1964,  p.  99-102. 

. J .  Taylor,  "Similarity  theory  in  the  relation  between  fluxes 
and  gradients  in  the  lower  atmosphere".  Royal  Met.  Soc.,  Quarterly 
Journal ,  Vol.  89,  1960,  p.  67-78. 

®G.E.  McVehil,  "Wind  and  temperature  profiles  near  the  ground 
in  stable  stratification".  Royal  Met.  Soc.,  Quarterly  Journal,  Vol.  90, 
1964,  p.  136-146. 
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have  failed  to  detect  any  systematic  variation  in  the  ratio  and  the 
question  has  remained  a  controversial  one. 

Air  Pollution  Measurements  in  Edmonton 

The  air  pollution  problem  in  Edmonton  is  not  yet  severe  but 

due  to  certain  climatic  factors  has  a  high  potential.  A.  Daniels 

in  1964  studied  the  problem  emphasizing  the  heat  island  and  its  rela- 

9 

tion  to  atmospheric  dispersion  .  In  his  study  he  concluded  that  there 
was  a  significant  increase  in  temperature  over  the  city’s  built-up 
areas  and  at  the  borders  of  the  central  business  district.  Over  the 
residential  areas  he  found  a  slow  increase  towards  the  city  center. 
Calculations  of  smoke  concentrations  based  on  a  modified  Sutton's 
diffusion  equation  supported  his  assumption  that  the  pollutants  were 
uniformly  distributed  vertically  throughout  the  heat  island.  The 
maximum  mixing  height,  computed  for  eight  winter  months , showed  a 
definite  variation  with  wind  direction,  being  highest  with  northerly 
winds,  lowest  with  southerly  winds.  This , together  with  the  fact  that 
southerly  winds  were  predominant  during  the  actual  months  of  his  in¬ 
vestigation,  convinced  him  that  any  plan  to  develop  polluting  industries 
to  the  south  of  the  downtown  area  should  be  ruled  out. 

In  the  same  year  the  Department  of  Health,  Government  of  Alberta, 
started  to  take  air  quality  measurements.  There  are  now  six  A.I.S.I. 
samplers  operating  on  two  hour  cycles  at  various  locations  in  the  city. 
The  pollution  data  have  been  correlated  with  wind  direction  and  wind 
speed  and  a  summary  of  the  results  has  been  published  on  a  yearly  basis. 

9 

A.  Daniels,  The  Urban  Heat  Island  and  Air  Pollution,  Unpublished 
M.Sc.  Thesis,  University  of  Alberta,  Edmonton,  1964. 
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The  most  favorable  meteorological  condition  for  the  accumulation 
of  pollutants  over  an  isolated  urban  area  would  be  a  combination  of 
persistent  light  winds  and  persistent  weak  vertical  mixing.  Two  para¬ 
meters  -  mixing  depth  and  mean  wind  speed  through  the  entire  air  volume  - 

have  been  proposed  as  basic  indicators  for  the  development  of  an  urban 

10 

pollution  potential  climatology  in  the  United  States  .  Extensive 
wind  profile  measurements  and  some  vertical  mixing  data  are  now  avail¬ 
able  for  many  rural  areas,  particularly  in  the  United  States.  Less 
well  known,  however,  are  the  effects  that  cities  may  have  on  mean 
winds  and  on  the  intensity  and  depth  of  vertical  mixing.  Diurnal 
changes  in  wind  speed  at  low  altitudes  can  be  used  as  a  relative  measure 
of  the  intensity  of  vertical  mixing  of  momentum  because  of  the  normal 
increase  of  wind  speed  with  height.  An  increase  in  v>ind  speeds  at 
low  altitudes  and  a  drop  of  speeds  at  higher  altitudes  results  in  an 
intensification  of  the  vertical  mixing.  An  investigation  done  by 
Geoscience  Research  Associates^  in  1967-68  for  the  city  of  Edmonton 
showed  that  while  the  average  wind  speeds  at  lOmeters  were  relatively 
strong  (8  to  11  m.p.h.  in  all  months),  periods  of  persistent  light 
winds  of  a  few  days  duration  were  not  uncommon,  particularly  in  winter . 
Such  periods  were  less  frequent  in  summer  because  of  intense  vertical 
mixing  that  results  in  a  pronounced  daytime  wind  speed  maximum.  Further- 


■^Report  of  Geoscience  Research  Associates  Limited  to  the 
Department  of  Health,  Government  of  Alberta,  Investigation  into  pollution 
in  Calgary  and  Edmonton,  April  1,  1967,  through  March  31,  1968. 


11 


Ibid. 
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more,  surface  and  upper- level  winds  in  Edmonton  provided  some  evidence 
of  urban  enhancement  of  vertical  mixing  in  midday  hours  in  winter  and 
at  all  hours  in  summer. 

Typical  weather  situations  in  which  the  soiling  index  would 
reach  relatively  high  levels  in  the  city  are  mainly  of  two  kinds: 

a.  Quasistationary  high  pressure  systems  centered  close  to  the  city, 

associated  with  low  winds  and  downward  air  movements  (subsidence). 

b.  Trapping  of  air  under  a  quasistationary  front  with  warm  air 

aloft  preventing  vertical  air  movement. 

The  Purpose  of  this  Study 

The  two  main  purposes  of  this  study  are  to  investigate: 

(1)  The  possibility  of  predicting  mean  daily  smoke  concentrations 
in  the  city  of  Edmonton; 

(2)  The  horizontal  and  vertical  temperature  distribution  over 
Edmonton,  and  the  vertical  temperature  distribution  over  Calgary. 

Method  of  Study 

Two  different  techniques  were  used  in  an  attempt  to  predict 
air  pollution  in  the  city  of  Edmonton.  To  arrive  at  useful  estimates 
of  the  mean  daily  smoke  output  in  the  city,  a  multiple  discriminant 
analysis  relating  smoke  output  with  some  meteorological  parameters  was 
applied  to  the  pollution  data  obtained  from  the  six  smoke  samplers. 

Using  regression  analysis  the  author  then  derived  and  tested  a  statistical 
model  for  the  variation  of  the  mean  daily  smoke  concentration  in  the  city. 

In  order  to  find  significant  features  of  the  horizontal  temper¬ 
ature  distribution  over  Edmonton  a  set  of  six  thermographs  were  placed 


xvii 


' 


■ 


in  two  concentric  circles  within  the  city.  The  data  thus  obtained  were 
then  analyzed  by  statistical  methods.  A  few  temperature  traverses  were 
undertaken  to  provide  more  detailed  data  on  selected  days.  A  temper¬ 
ature  traverse  is  a  series  of  temperature  measurements  taken  along  a 
pre-determined  route.  In  an  attempt  to  find  significant  features  of 
the  vertical  temperature  distribution  in  the  cities  of  Edmonton  and 
Calgary  thermometers  at  different  altitudes  were  mounted  on  towers. 

The  data  were  analyzed  statistically  and  the  results  were  compared  with 
results  from  other  sites. 
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CHAPTER  I 


INSTRUMENTATION 


General 

This  chapter  describes  the  various  instruments  used  in  the 
author's  investigation.  As  data  are  being  collected  from  both  Edmonton 
and  Calgary,  the  instruments  from  both  sites  are  described.  However 
the  analysis  will  be  focused  on  the  Edmonton  data. 

The  CN  Tower  and  Calgary  Tower  Instruments 

The  vertical  temperature  structure  in  Edmonton  is  measured  by 
electrical  resistance  thermometers^.  The  system  consists  of  a  balanced 
bridge,  recorder  and  a  sensing  element,  functioning  as  follows: 

The  Circuit. 

The  measuring  element  forms  one  leg  of  the  bridge.  When  the 
resistance  of  the  measuring  element  changes,  current  flows  in  the 
bridge  and  a  voltage  unbalance  is  created  in  the  circuit.  This  voltage 
is  sufficiently  amplified  to  operate  a  balancing  motor.  The  balancing 
motor  adjusts  the  sliding  contact  on  a  slidewire  to  restore  electrical 
balance  in  the  bridge  circuit  and  simultaneously  moves  the  pen  to 
record  on  the  chart  the  variable  being  measured.  The  balancing  action 
starts  with  the  slightest  change  in  the  resistance  and  so  the  measuring 
procedure  is  practically  instantaneous. 

Sensing  Element. 

The  system  is  equipped  with  a  resistance  thermometer  bulb,  a 
temperature  sensing  element  which  changes  resistance  with  variations 


■^Bristol  Model  66A  12064,  Waterbury  20,  Conn. 
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in  temperature.  Heating  of  the  extension  wire  is  supposedly  compensated 
for.  Thus  the  accuracy  of  temperature  measurements  at  the  bulb  are 
not  affected  by  temperature  differences  or  changes  anywhere  along  the 
wires.  Resistance  thermometers  are  widely  used  because  of  the  follow¬ 
ing  advantages;  remote  recording,  high  sensitivity,  small  dimensions 
and  fast  response.  The  thermometers  are  metallic  and  usually  made  from 
platinum,  nickel  or  copper.  The  relation  between  resistance  and  temper¬ 
ature  for  the  bulb  can  be  expressed  as : 


R  =  Ro  (1  +  aT  +  bT 


2 


+ - ) 


(1) 


Where  Ro  is  the  resistance  at  0°C,  T  is  temperature  in  °C,  and  a  and  b 
are  constants.  This  relation  is  almost  linear  if  the  temperature 
is  small.  If  Eq .  (1)  is  replaced  by  the  linear  relation: 


R  =  Ro  (1  +  aT) 


(2) 


it  can  be  shown  that  the  maximum  temperature  error  in  the  temperature 

2 

interval  10-40°C  is  less  than  0.4°C.  This  value  applies  strictly  to 
a  platinum  bulb  but  a  similar  one  can  be  derived  for  any  other  metal. 

3 

Electrical  resistance  thermometers,  manufactured  by  the  Foxboro 
Company  were  mounted  on  the  Calgary  Tower.  The  sensors  are  housed  in 
aspirated  shields  for  adequate  ventilation  and  protection  from  direct 
sunlight.  The  recorder  is  of  a  balanced  bridge  type. 


o 

W.E.K.  Middleton  and  A.F.  Spilhaus,  Meteorological  Instruments, 
University  of  Toronto  Press,  Toronto  1953,  p.  84. 

q 

Model.  No.  DB-23S-26W,  Foxboro  Co.,  Lasalle,  Que. 
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The  A.I.S.I.  Smoke  Sampler 

The  equipment  used  is  an  automatic  spot  sampler  developed  by 
the  American  Iron  and  Steel  Institute  .  In  this  instrument  air  is  drawn 
through  a  filter  tape  which  is  held  firmly  by  a  clamp.  After  the  re¬ 
quired  sampling  time,  in  this  case  two  hours,  the  clamp  is  released 
and  a  fresh  area  of  filter  is  moved  into  the  sampling  head.  The  soiled 
spots  are  evaluated  by  measuring  the  optical  density  (O.D.)  of  the 
spot  with  a  photometer.  This  is  set  to  read  1007,  transmittance  through 
a  clean  filter.  A  value  can  then  be  read  for  the  soiled  spot.  The 
optical  density  can  be  written: 

O.D.  =  log  Io  (3) 

I 

Where  Io  is  the  intensity  of  light  transmitted  through  clean  tape  and 
I  is  the  intensity  of  light  through  soiled  tape.  As  this  is  a  non¬ 
linear  relation  caution  must  be  observed  in  interpreting  samples  which 
fall  below  about  507o  light  transmission.  Above  this  value  the  relation 
is  approximately  linear.  For  comparison  with  samples  taken  from  other 
places  the  optical  density  is  converted  to  a  unit  called  the  COH  unit 
defined  as: 

1  COH  unit  =  100  x  O.D.  (4) 

COH  is  an  abbreviation  for  "coefficient  of  haze".  The  density  of  the 
filter  stain  depends  on  the  volume  of  air  drawn  through  the  tape. 
Therefore  the  volume  is  divided  by  the  area  of  the  sampled  spot  to 

4 

Model  No.  F2 ,  Research  Appliance  Co.,  Allison  Park,  Penn. 
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obtain  a  linear  unit  of  air. 

L  =  t  x  Q  (5) 

A 

Where  t  is  the  sampling  time  in  minutes, 

Q  is  the  average  flow  rate  for  the  sampling  time  in  (feet)  min“^, 
and 

2 

A  is  the  sampling  area  in  (feet)  . 

For  convenience  the  quantity  of  air  is  expressed  in  thousands  of 
linear  feet. 

L  _  _ t  x  Q  (6) 

1000  x  A 

Thus : 

COH  units  per  1000  linear  feet_  100  x  O.D.  _  10~*  x  O.D.  x  A  (7) 

L  t  x  Q 

For  the  samplers  used  in  this  investigation  the  following  values 
are  valid: 

t  =  120  minutes 
Q  =  0.25  cfm 

2 

A  =  0.0054  ft.  (1  inch  diameter) 

Substituting  these  values  in  Eq.  (7)  gives: 

COH/lOOO  ft  _  105  x  O.D.  x  0.0054  =  18  x  O.D.  (8) 

120  x  0.25 
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Thermographs 

Seven  thermographs  were  installed  in  Edmonton  and  Calgary  for 
measuring  the  horizontal  temperature  distributions.  These  instruments 
employ  a  bimetal  element  consisting  of  a  compound  strip  of  metal,  form¬ 
ed  by  welding  together  two  bars  of  different  metals.  The  resulting 
strip  bends  as  its  temperature  changes,  because  of  the  difference  in 
expansion  of  the  two  metals.  The  deformation  is  then  transmitted  through 
a  linkage  system  to  a  recording  pen.  The  data  are  recorded  on  a  seven 
day  clock  and  drum  assembly  whose  chart  is  printed  in  one  degree  intervals 
over  a  range  of  -40°C  -  +40°C.  The  instrument  cover  which  is  finished 
in  white  enamel  provides  shielding  from  solar  radiation  while  allowing 
free  movement  of  air  around  the  sensing  element.  Before  being  put  into 
operation  all  instruments  were  carefully  calibrated  against  an  aspirated 
Assman  psychrometer  and  yielded  a  maximum  error  of  ±0.3°C  in  the  temper¬ 
ature  range  of  -26°C  -  +19°C.  The  instruments  are  at  present  possibly 
reading  up  to  6-8°C  too  high  under  conditions  of  light  wind  and  bright 
sunshine.  These  large  errors  are  attributed  to  inadequate  ventilation. 
Other  large  errors  occasionally  resulted  from  acts  of  vandalism  and 
accidental  damage  to  the  instruments.  In  the  absence  of  strong  solar 
radiation  comparisons  with  standard  instruments  showed  differences  of 
±0.5°C  or  less. 

Mobile  Unit 

Instrument  used  for  the  Temperature  Traverses 

In  addition  to  the  thermographs,  temperature  traverses  were 
taken  to  obtain  a  more  complete  picture  of  the  horizontal  distribution 

Thermograph  No.  0T04,  Eiwa  Seiko  Company,  Osaka,  Japan, 
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of  temperature.  As  the  mobile  unit  is  moving, it  was  necessary  to 
have  an  instrument  that  responded  quickly  to  small  changes  in  temperature. 
It  was  thought  that  an  electrical  resistance  thermometer  would  meet 
this  requirement.  Therefore  such  a  thermometer,  manufactured  by  the 
same  company  as  the  Calgary  Tower  thermometers,  was  mounted  on  top  of 
a  van  belonging  to  Geoscience  Research  Associates  Limited.  The  instru¬ 
ment  was  shielded  from  direct,  sunlight.  Sufficient  ventilation  is  pro¬ 
vided  for  when  the  van  is  moving.  The  sensor  is  connected  with  a 
recorder  placed  inside  the  van  where  it  can  be  easily  read  by  the 
driver . 

The  Lag  Coefficient  of  the  Thermometer 

When  a  thermometer  of  any  sort  at  a  certain  temperature  is 
placed  in  a  medium  at  a  different  temperature  it  does  not  immediately 
take  up  the  temperature  of  the  medium  but  approaches  it  asymptotically 
at  a  rate  which  depends  on  the  materials  and  dimensions  of  the  thermometer 
and  also  on  the  properties  of  the  medium  itself. 
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The  change  in  the  indication  of  the  thermometer  with  time  is 
given  by  the  differential  equation: 


di 

dt 


(9) 


where  i  is  the  indicated  value,  r  is  the  true  value,  and  a  is  the  lag 
coefficient  with  dimensions  of  time. 

Let  us  now  assume  that  a  is  a  real  constant  and  also  that  the 
temperature  of  the  medium  does  not  change.  Then  integration  of  Eq . 

(9)  gives: 

fdt  <10> 

-  r  dt  a 

'  o 

where  i  is  the  value  of  i  at  time  zero 
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(13) 


If  the  substitution  t= a  is  done  in  Eq.  (13)  we  see  that  the 
lag  coefficient  is  the  number  of  seconds  required  for  the  difference  of 
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temperature  to  be  reduced  to  e"1  =  l/e  =  1/2.718  =  37%  of  its  initial 
value . 

The  lag  coefficient  a  is  highly  dependent  on  ventilation  and 
decreases  rapidly  with  increasing  ventilation.  As  the  lag  coefficient 
is  not  known  in  this  case  it  has  to  be  calculated.  According  to  Newton's 
law  of  cooling  the  rate  of  change  of  temperature  (0)  for  the  thermometer 
can  be  written: 

K(0-T)  =  -MC  (14) 

where  M  is  the  mass  of  thermometer,  C  is  the  specific  heat  of  the 
thermometer,  and  K  is  the  coefficient  of  dissipation,  and 

T  is  the  air  temperature 

Eq .  (14)  and  Eq .  (9)  combined  gives: 


_d9  _  K  (Q-T) 

dt  MC 


d0  =  _  1 

dt  a 


(0-T) 


Then : 

a  =  M2 
KL 


(15) 


where  L  is  the  length  of  the  probe 

Thus  the  lag  coefficient  contains  the  dissipation  coefficient  of 
the  thermometer  as  a  factor  and  is  also  a  function  of  the  specific  heat 
of  its  material. 

% 

The  bulb  is  made  of  brass  coated  with  nickel  which  is  estimated 
to  be  not  over  3%  of  the  total  weight.  The  bulb  is  cylindrical  with 
diameter  0.635  cm.  and  length  26  cm.  The  difference  in  density  for  the 

_  o 

two  metals  is  less  than  0.4  gem  ,  the  difference  in  specific  heat  is 
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0.018  cal.  g.“^-  deg.”-*-  and  the  difference  in  the  dissipation  coeffi¬ 
cient  is  0.16  cal.  sec.’l  deg.~^  cm.  ^ .  In  spite  of  the  large  differ¬ 
ence  in  the  dissipation  coefficient  the  author  chose  to  base  all  cal¬ 
culations  on  brass  alone.  The  mass  of  the  bulb  can  be  calculated  from: 


M  =  PV  =  PttR2L  (16) 

which  gives  M  =  70  g. 

The  specific  heat  for  brass  at  20°C  is  0.09  cal.  g.-^,  its  density 
P  “8.5  g.cm.'^and  its  dissipation  coefficient  at  the  same  temperature 
is  0.30  cal.  sec.’-*-  deg.“^-  cm.'^. 

The  lag  coefficient  a  can  then  be  calculated  according  to 
Eq.  (15): 


70  x  0.09 
0.3  x  26 


0.7  sec. 


Thus  the  thermometer  will  register  63%.  or  l/e  of  a  sudden 
temperature  change  in  0.7  sec. 

As  an  example  of  the  thermometer's  efficiency  consider  a  sudden 
temperature  change  of  1°C.  How  far  can  the  van  run  before  90%  of  the 
change  is  recorded? 

The  time  it  would  take  to  do  this  can  be  obtained  from  Eq.  (13): 

=  JL 

10 

e1/"  =  10 

and 

t  =  a In  10  =  2.3  a  =  1.61  sec. 
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If  the  van  is  driven  at  maximum  speed  30  m.p.h.  the  distance 
covered  would  be: 

30  x  1.609  x  103  x  1.61  ,  22  m  (72  ft.) 

36  x  102 

This  maximum  distance  is  assumed  to  be  small  enough  for  the 
traverses . 

Location  and  Physiography  of  the  Edmonton  Region 

Edmonton  with  a  population  of  394,000  (1968)  and  an  area  of 
85  square  miles  (1966)  is  the  capital  of  Alberta,  located  at  latitude 
53°33’  North  and  longitude  113o30I  West.  The  valley  of  the  North  Sask¬ 
atchewan  is  the  main  channel  in  the  area.  It  cuts  across  the  city  from 
southwest  to  northeast.  The  river  valley  is  deep  and  twisting,  ranging 
in  width  from  half  a  mile  in  the  southwest  to  a  mile  in  the  center  of 
the  city.  The  average  depth  is  160  feet,  but  in  places  it  is  over  200 
feet  deep. 

The  angle  of  slope  of  the  river  valley  varies  from  5°  in  the 
university  area  to  45°  at  Grierson  Hill. 

Several  tributary  ravines  enter  the  main  valley  from  both  sides. 

The  residental  area  of  the  city  is  mainly  of  low  density 
(  <50  persons  per  net  acre)  with  exception  of  the  city  core  which  can 
be  classified  as  being  of  medium  density  (50-89  persons  per  net  acre)  . 
The  houses  in  these  areas  are  generally  one-  and  two- storey  homes  and 
are  rather  close  to  one  another.  The  Commercial  Districts  consisting 

6From  City  Land  Use  Map,  1965,  City  Planning  Office,  Edmonton. 
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of  closely  spaced  and  taller  buildings  are  located  in  the  core  of  the 
city,  and  also  along  82nd  Avenue  and  in  Jasper  Place.  Industries  are  con¬ 
centrated  to  the  east,  south  and  north-west.  There  are  several  open  spaces, 
such  as  ravines,  parks,  recreation  areas,  sports  fields,  playgrounds 
and  an  airport,  amounting  to  roughly  12%  of  the  total  built  up  area, 

(Fig.  1). 

Location  of  the  Edmonton  Instruments 
CN  -  Tower 

Electrical  resistance  thermometers  at  three  different  levels 
(  15  m. ,  57  m.  and  113  m.)  were  installed  on  a  large  office  building 
(CN  Tower)  located  in  downtown  Edmonton.  The  lowest  and  highest  were 
mounted  on  8  ft.  booms  projected  from  the  building  on  the  north  wall. 

The  sensors  were  housed  in  aspirated  shields  to  provide  good  ventilation 
and  protection  from  direct  solar  radiation.  The  intermediate  level 
thermometer  was  placed  in  an  air  intake  on  the  east  wall  of  the  build¬ 
ing.  Data  collection  was  started  on  a  continuous  basis  on  June  1,  1967. 
Calibration  carried  out  at  monthly  intervals  showed  no  drift  or  loss 
of  accuracy.  However,  a  study  of  vertical  temperature  profiles  for 

July  and  December,  1967,  showed  some  evidence  to  indicate  that  the 

7 

observed  values  at  57  m.  were  not  representative  . 

Thermographs 

The  locations  of  the  seven  thermographs  had  to  fulfill  at 
least  the  following  criteria: 

a.  They  had  to  be  placed  on  open  areas  away  from  buildings  to 

^Geoscience  Research  Associates  Limited,  Annual  Report : 
Investigation  into  Air  Pollution  in  Calgary  and  Edmonton.  Edmonton, 

1968^ 
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minimize  the  effect  of  heat  radiation. 

b.  At  the  same  time  they  had  to  be  placed  sufficiently  close  to  a 
continuously  attended  building  to  minimize  damage  by  vandalism. 

c.  Last  but  not  least,  they  had  to  be  distributed  in  such  a  way 
as  to  effectively  measure  the  temperature  differences  as  one 
approaches  the  city  center  from  the  outskirts  and  so  give  an 
indication  of  the  intensity  of  the  heat  island. 

The  following  places  were  selected:  (Fig.  2) 

#1  1160  River  Valley  Road 

#2  City  Hall,  100  Street  and  103  Avenue 

#3  University  of  Alberta  Campus 

#4  11504  -  37 B  Avenue 

#5  Jasper  Place  High  School 

#6  12735  -  101  Street 

#7  6625  -  101  Avenue 

Continuous  temperature  records  were  started  in  January,  1968. 
As  unit  #3  was  often  subjected  to  vandalism  it  was  moved  to  11333  -73 
Avenue  on  October  3,  1968. 

Air  Samplers 

Since  August,  1964,  there  have  been  six  A.I.S.I.  samplers 
operating  in  the  city,  all  of  which  were  installed  by  the  Department 
of  Health  of  the  Province  of  Alberta.  The  locations  of  these  stations 
are:  (Fig  3) 

#1  CN  Tower 

#2  Administration  Building  (commercial) 

#3  117  Avenue  -  147  Street  (industrial) 
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Location  of  Edmonton  Thermograph  Stations 


dufO  o>ofl*l( 

rntq  «,jC  ;oL 

-15- 


00 «  x. 


\  Sherwood 

r:"‘ 

>Park 


AIR  POLLUTION 


DATA 


COLLECTING 


STATIONS 


MIL  ES 


Figure  3 

Edmonton  Air  Pollution  Data  Collecting  Stations 
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#4  Sherwood  School  (residential) 

#5  C.N.R.  North  Edmonton  (industrial) 

#6  97  Street  80  Avenue 

A  data  summary  is  published  each  month  by  the  Department  of 

Health. 

Location  of  the  Calgary  Instruments 

Tower 

The  data  tower  in  Calgary  is  located  at  the  west  end  of  Bonny- 
brook  Bridge  in  an  open  area  which  is  grass  ^covered >at  an  elevation  of 
3,345  feet  above  sea  level.  (Fig.  4)  The  tower  is  of  three-cornered 
tubular  metal  construction  measuring  18  inches  on  each  side  and  300 
feet  high.  Electrical  resistance  thermometers  were  installed  at  four 
different  levels  (10  ft.,  100  ft.,  200  ft.,  and  300  ft.).  As  on  the 
CN  -  Tower,  each  sensor  is  housed  in  an  aspirated  shield.  Continuous 
data  collection  was  initiated  on  March  1,  1968. 

In  the  initial  stage  the  300  feet  sensor  showed  considerable 
drift,  due  to  poor  ventilation.  This  was  corrected  on  September  4, 


1968. 
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Figure  4 


Location  of  Calgary  Instrument  Tower 


CHAPTER  II 


PREDICTION  OF  SMOKE  IN  EDMONTON 


General 

This  chapter  deals  with  the  problem  of  predicting  pollution. 

Two  different  procedures  were  used  to  obtain  prediction  equations 
for  a)  smoke  emission  rates  (Q)  and  b)  absolute  concentrations 

(  X). 

a)  Prediction  of  Smoke  Emission  Rates  (Q) 

Although  the  author  was  aware  of  the  difficulties  in  obtain¬ 
ing  useful  estimates  of  smoke  emission  rates  using  meteorological 
variables,  it  was  nevertheless  decided  to  apply  a  multiple  discriminant  an¬ 
alysis  to  the  data.  Similar  work  has  been  done  In  the  United  States  by  e.g. 
M.E .  Miller  and  G.C.  Holzworth*"  who  derived  equations  for  source  strength 
versus  mean  daily  temperatures  for  three  American  cities.  The  results 
they  obtained  in  testing  the  formulae  were  very  encouraging.  In  all 
cases  the  simple  correlation  coefficient  between  observed  and  calcu¬ 
lated  values  was  higher  than  0.8. 

In  this  study  the  derived  discriminant  function  failed  to 
predict  the  output  of  smoke,  the  reason  being  the  choice  of  predictors. 
Temperature  variations  did  not  account  for  a  significant  portion  of 
the  variation  of  smoke  output  in  Edmonton.  A  possible  explanation  for 
this  is  that  space  heating,  which  is  temperature  dependent,  is  not  a 
major  source  of  smoke,  due  to  the  use  of  natural  gas.  The  slight 


■^M.E.  Miller  and  G.C.  Holzworth,  An  Atmospheric  Diffusion 
Model  for  Metropolitan  Areas.  Air  Resources  Field  Research  Office, 
ESSA  Cincinnati ,  Ohio,  1967. 
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variation  with  temperature  in  output  from  combustion  engines  is  not 
enough  to  be  significant.  A  description  of  the  analysis  that  leads 
to  the  three  prediction  equations  is  given  below. 

The  technique  of  discriminant  analysis  was  originated  by  R.A. 
Fisher  in  1936  and  later  extended  by  J.G.  Bryan  to  multiple  variables. 
The  method  can  be  described  as  a  multivariate  statistical  model  which 
expresses  a  set  of  mutually  uncorrelated  discriminant  functions  as 
linear  functions  of  a  fixed  set  of  predictors  (independent  variables). 
The  functional  relationships  are  determined  from  a  sample  such  that 
the  ratio  of  the  between-group  sum  of  squares  to  the  within-group  sum 
of  squares  is  maximized  for  each  discriminant  function.  In  other 
words  a  criterion  that  maximizes  the  distance  between  the  group  means, 
and  minimizes  the  spread  of  the  points  about  the  means,  i.e.  a 
criterion  that  maximizes  the  ratio: 


SSB  (Y  j  ) 
SSW(Yj) 


where 


SSB  (Y  j  ) 


and 


SSW  (Y  . ) 
J 


£ 

k=l 


(1) 


(2) 


(3) 


k=l, - n, 

’  g 


g=l,  ~~ G 
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Mat  hematic  a  1  Procedure 


Predictor  Space 


Let  there  be  G  mutually  exclusive  and  exhaustive  groups  defined 
for  the  predictand  variable.  Assume  there  are  P  predictors  denoted  by 

Xp  (p  =1,  ---  P). 

Let  the  number  of  observations  in  each  group  of  the  dependent 
sample  be  n  (g  =  1 ,  - G) 


Let 


N  _  E 


n 


8=1 


g 


(4) 


Determine  the  quantities: 

The  pooled  within-group  sum  of  squared  deviations  about  the 


group  means, 


SSW(Xp) 


G 

E 

8=1 


ng 

£ 

k=l 


(Xpgk  -  xpg) 


(5) 


(p  =  1,  P) 


and  the  sum  of  squared  deviations  between  group  means  and  the  grand  mean, 


SSB(Xp) 


G 

E 

g=l 


n 


(X 


Pg 


Xp) 


(p  =  1,  ---  p) 


(6) 


Also  required  are  the  sum  of  products  within  (SPW)  and  between  (SPB)  groups 
For  example,  for  predictors  Xp  and  X^ 


(P,q  =  1,  P  5  where  p  j  g) 

r  nrr 

G  g 

SPW  (X  X  )  =  E  E  (X  k  -  Xpg)  (X  k  -  Xqg) 

* 2 *  4  g=l  k=l 

(p >q  =  i,  ---  p;p  f  q) 


(7) 


2This  section  is  a  summary  of  R.G.  Miller's  article,  "Statistical 

Prediction  by  Discriminant  Analysis,"  Meteor.  Monographs,  VoL.  4,.Amer. 

Meteor.  Soc. ,  1962. 
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and 

SPB (XpXq)  = 


J.  VV)  -  V  (Xpq  '  V 


(p,q  =  l,  —  p  ;  p  f  q) 


(8) 


A  pooled  within-group  matrix  W  is  constructed  from  the  derived 


quantities 


W  = 


< 


f  ssw(xx)  spw(XjX2)  ---  spwcx^p.x)  spwo^Xp) 

SPW(X1X2)  SSW f(X2)  ---  SPW (X2Xp_i)  SPW(X2X  ) 


SPW (X][Xp_ ]_)  SPW(X2Xp_1)  --  SSW(Xp„1)  SPW(Xp_1Xp) 
v  SPWCX^p)  SPW(X2Xp)  ---  SPW(Xp_1Xp)  SSW(Xp) 


/ 


(9) 


B  = 


Similarly  a  between-group  matrix  B  is  constructed  as  follows: 


< 


/SSB(X1)  SPB(X1X2) 

spb(x1x2)  ssb(x2) 


SPB  (XLXp_  x)  SPB(XxXp) 
SPB(X2Xp_!)  SPB(X2Xp) 


\ 


SPB(X1Xp.1)  SPB(X2Xp.1)  ---  SSBCXp.p  SPB(Xp_LXp) 
VSPBCXpXp)  SPB(X2Xp)  —  SPB(Xp_]X  )  SSBCXp)  J 


(10) 


Discriminant  Space 

The  multiple  discriminant  analysis  procedure  makes  use  of 

both  W  and  B  to  determine  a  new  set  of  functions  Yj  {  j  =  1, - min. 

(G  -  1,  P)  },  which  are  linear  in  the  original  predictors  X  .  The 
method  involves  1)  an  inversion  of  W,  2)  a  pre-multiplication  of  B 
by  the  inverse  of  W,  and  3)  a  determination  of  the  eigenvalues  and 
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characteristic  vectors  of  the  resulting  matrix  W~^B. 

This  is  obtained  by  equating  to  zero  the  derivative  of  Eq.  (1) 
which  is  to  be  maximized. 


Thus 


S SW (Y  .  )  (SSB  (Y  ; )  )  -  SSB  (Y:)  (SSW(Y,)) 

J  3Vjp  J  J  3  V._  J 


JP 


<\p 

3  V.  J 
JP 


{  SSW (Y j )  } 


=  0  (11) 


By  definition: 


Y.  .  =  Z  V.  XncrV 

jpk  JP  P§k 


(12) 


Differentiating  Eq .  (2)  with  respect  to  V ^  gives: 


i_  (SSB  (Y  j) )  =  2  E  ng  (Yjg.-Yj.) 


3  Y  . 
Jg 


3vjp 


g 


3V 


JP 


8Yj.. 


(13) 


where  one  dot  means  summing  over  k  and  two  dots  mean  summing  over  k 


and  g. 


From  Eq.  (12) 


9YJP 


3  V . 
JP 


"  ?  XPg  -  V 


(14) 


9  Y  • 

- — '  =  z  Z  X  =  X 

3V.  k  g  P  p' 

JP 


Substituting  Eq.  (14)  into  Eq .  (13),  again  making  use  of  Eq. 


(12)  gives: 


-23- 


L_  (SSB  (Y  j ) )  =2  E  n  (Z  V.  X 


3  V . 
JP 


iT~0  e  v.  x  }  (x  -  L  ) 

g  -  p  JP  pg-  JP  p. ♦  pg.  p- • 


2  E  n  E  V.  (X  _X  )  (X  -X  ) 
g  jp  v  pg.  p. . ;  v  qg-  q- 
g  p 


(15) 


where  q  is  a  dummy  index. 

This  last  expression  is  by  definition  equivalent  to  (B)  V 


JP 


Similarly : 


3_  (SSW(Yj)) 


3V. 

JP 


(W)  V. 

V  '  JP 


2  E  E  E  V-  (X  v)-  (X  )  (X  ,  -  X  )  = 

JPV  Pgk/  v  Pg*/V  qgk  qg.' 

g  k  p 


(16) 


Eq.  (11)  can  now  be  written 


3  SSB  (Y  . ) 

3  (X.)  =  SSW  (Y  •)  _ -  SSW  (Y  •)  X. 

J  J  5  T7  J  J 


3  V . 
jp 


3  V  • 
JP 


3  SSW (Y . ) 


3  V . 
JP 


=  0  (17) 


Thus  : 

3  SSB (Yj ) 

3Vjp 


3  SSW (Y j ) 


3  V . 
JP 


=  0 


(18) 


which  with  Eq .  (16)  and  Eq .  (18)  gives 


(B  -  XjW  )  Vjp  =  0 


(19) 


The  eigenvalues  (X  )  are  then  obtained  by  solving  the  deter- 


minantal  equation 


W"1  B  -  X  I  I  =  0 


(20) 


' 


■ 
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where  I  is  the  unit  matrix. 

The  characteristic  vectors  (Vj)  which  represent  the  weights 
on  the  predictors  are  solutions  of  the  equations: 

{W-1  B  -  Xi  I  }V-  =  0 

J  J  (21) 

(j  =  1,  -  min(G  -  1,P)) 

where  I  is  the  unit  matrix. 

The  discriminant  functions  can  now  be  written  as 

Y .  =  Vu  Xi  +  X?  + - V.  X 

J  J1  1  Jz  1  JP  P 

(22) 

(j  =  1,  -  min(G  -  1,P)) 

Predictors 

In  this  study  three  predictors  were  selected,  the  selection 
being  based  on  the  ready  availability  of  the  data  and  on  some  experience. 
All  three  predictors  were  transformed  into  dummy  variables  to  simplify 
calculations.  Below  follows  a  description  of  the  predictors  and  the 
values  they  take  when  transformed. 

Predictor  X-^ :  Above  or  below  monthly  normal  minimum  temperature  in 
°F.  (based  on  a  thirty  year  average) .  X^  is  assigned  a  1  if  less  than 
or  equal  to  the  monthly  normal  minimum  temperature,  and  0  if  it  is 
greater  than  this  temperature. 

Predictor  X^  :  This  predictor  also  makes  use  of  monthly  normal 
minimum  temperature.  is  assigned  a  1  if  it  is  either  much  below  or 

much  above  monthly  normal  minimum  temperature,  and  a  0  if  it  is  either 
slightly  below  or  slightly  above.  Much  above  and  much  below  were 
defined  as  normal  minimum  temperature  plus  10° F  or  more  and  normal  min¬ 
imum  temperature  minus  10°F  or  more,  respectively. 


■ 

. 
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Predictor  X^:  Wind  direction.  was  assigned  a  1  if  the 

wind  direction  was  from  NW  -  N  -  NE ,  and  a  0  in  all  other  cases. 

Daily  values  were  taken  for  six  winter  months,  October  through  March 
(1967-68)  and  compared  with  the  predictand.  The  data  were  readily 
available  from  the  Monthly  Meteorological  Summary  prepared  by  the 
Department  of  Transport  (Canada)  for  the  city  of  Edmonton. 

Predictand 

The  predictand  was  chosen  to  be  smoke  emission  rate,  with 

-2  -1 

units  g.  cm.  sec. 

As  smoke  emission  data  were  not  available  they  had  to  be 

calculated  from  a  mathematical  model  relating  smoke  emission  rate  with 

concentration  of  smoke,  mixing  height,  and  wind  speed  as  follows: 

Assuming  a  uniform  smoke  concentration  within  a  box- like  volume, 

and  assuming  a  constant  smoke  production  rate;  the  following  result  can 
3 

be  derived: 

X  _  j-  X  0  (hp)  -  QL  -i  exp  (“lit )  QL  (23) 

h  uh  L  uh 

When  t  approaches  large  values  the  concentration  will  reach  a 
state  of  equilibrium: 

X  =  QL  (24) 

uh 

Solving  Eq.  (24)  for  Q  gives: 

Q  =  uh  X  (25) 

L 

3S.  Djurfors,  "  An  urban  pollution  model".  The  Albertan  Geographer, 
No.  4,  University  of  Alberta,  1968  P.  65-68. 
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where  u  is  the  mean  wind  speed  in  cm.  sec.-^  in  the  layer  under 
consideration,  h  is  the  mixing  height  in  cm.,  y  is  the  ground  con¬ 
centration  of  smoke  in  g.  cm.'^,  and  L  is  the  effective  cross-city 
dimension  parallel  to  the  wind  direction  and  varying  with  the  wind 

speed.  In  this  study  L  was  taken  to  be  constant  with  the  value  128 
4 

x  10  cm.  Assuming  the  city  lapse  rate  to  be  dry  adiabatic  due  to 
vigorous  turbulent  mixing,  the  mixing  height  is  obtained  as  the 
intersection  between  the  dry  adiabat  from  the  observed  ground 
temperature  in  the  city  and  the  observed  temperature  profile  over 
the  country^.  This  can  be  stated  mathematically  as  follows: 

Let  the  country  lapse  rate  be: 


dT  = 
dZ 


-  Y 


(26) 


and  let  the  city  lapse  rate  (T  —  dry  adiabatic)  be: 


dZ 

Integration  from  Z=0  to  Z=h  and  substitution  gives 
Tx  -  To  =  h  (r^>  =  ha 


(27) 


(28) 


Thus 


h  = 


T  -  T 
Xx  o 


a 


(29) 


where  Tx  is  the  city  ground  temperature  and  TQ  is  the  country  ground 


^P.W.  Summers,  An  Urban  Ventilation  Model  Applied  to  Montreal. 
Ph.  D.  thesis,  Montreal,  1964,  P.  106-109. 
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temperature.  The  concentration  of  smoke  (X  )  is  given  in  COH-units 
which  is  a  relative  measure  and  so  has  to  be  transformed  into  metric 
units.  Sullivan  found  that  the  best  expression  for  the  calculation 
of  smoke  density  from  smoke  stains  of  one- inch  diameter  was  as  follows: 


where  M 

-3 

yg.rn  . 
equal  to 


M1  =  a  (D  x  100)b 

0.67 

V 


(30) 


is  mass  of  particles  with  diameter  less  than  10  microns  in 
Here,  D  is  the  optical  density  by  reflection,  a  is  a  constant 
33.0,  V  is  air  volume,  and  b  is  another  constant  equal  to 


1.24. 


Thus : 


M 


33.0  (D  x  100)1*24 
(120  x  0.25)0,67 


3.38  (D  x  100)1*24 


(31) 


The  optical  density  by  reflection  (D)  is,  however,  related  to 
the  optical  density  by  transmission  (O.D.)  in  the  following  way: 


D  =  43  x  O.D.  =  0.98  x  O.D.  (32) 

44 

From  Chapter  I  (Eq.  8): 

O.D.  =  (  COH/1000  ft  0/18.0 


Substitution  of  the  above  relation  and  Eq .  (32)  into  Eq .  (31) 


gives : 


M1  =  3.38(98/18.0  x  COH/1000  ft.)1,24  = 

-I  r\  r 

3.38(5.44  x  COH/1000  ft.)  =27.6  (COH/1000  ft.)1*24 


M  refers  only  to  the  fraction  of  particles  with  diameter  less  than 


10  microns.  To  involve  all  particles  it  is  assumed  that  the  fraction 


•*J.L.  Sullivan,  The  Calibration  of  Smoke  Density.  Presented  at 
the  55th  annual  meeting  of  Air  Pollution  Control  Association,  Chicago, 
1962. 
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of  small  particles  on  the  average  is  707.  of  the  total  mass  of  particles. 

ThuS:  1  9 A 

M  =  39.5  (COH/IOOO  ft.)  (34) 

where  M  is  the  total  mass  of  particles  in  yg.  m-3. 

Hourly  readings  of  smoke  concentrations  were  obtained  from  the 
Department  of  Health  for  each  of  the  six  sampling  stations  and  were 
averaged  over  24  hours.  This  was  done  to  smooth  the  data  as  previous 
tests  had  shown  that  2-hour  and  even  8-hour  averages  were  poorly  cor¬ 
related  from  one  station  to  another.  In  Table  1  are  shown  the  simple 
correlation-coefficients  with  each  pair  of  stations  for  24  hour  read¬ 
ings  averaged  over  six  winter  months,  October  through  March,  1967-68. 


Table  1 

Correlation  Matrix  of  Smoke  Sampling  Stations 


Stn.  No. 

Correlation 

Coefficients 

1 

2 

3 

4 

5 

6 

1 

X 

0.44 

0.29 

0.17 

0.23 

0.19 

2 

0.44 

X 

0.48 

0.45 

0.45 

0.47 

3 

0.29 

0.48 

X 

0.29 

0.62 

0.30 

4 

0.17 

0.45 

0.29 

X 

0.33 

0.41 

5 

0.23 

0.45 

0.62 

0.33 

X 

0.41 

6 

0.19 

0.47 

0.30 

0.41 

0.41 

X 

The  values  are  still  quite  small  and  it  may  be  better  to  use 
longer  time  averages  in  the  future.  In  spite  of  all  this  it  was  de¬ 
cided  to  use  the  24-hour  readings.  These  readings  were  then  averaged 
over  all  stations  to  obtain  a  value  supposedly  representing  the  average 
daily  smoke  concentration  in  Edmonton.  Making  use  of  Eq .  (25)  the 

mean  daily  smoke  emission  rate  (Q)  was  calculated.  The  predictand  was 


•  fV 

.  > 
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divided  into  three  classes  as  follows: 
Q3  =  (0  -  0.2)  x  10'10 


g. 

cm . 

-2 

sec . 

-1 

-10 

-2 

-1 

g. 

cm. 

sec . 

-10 

-2 

-1 

g- 

cm. 

sec . 

Q3  =  (0.501  -  6.0)  x  10 

Discriminant  Functions 

The  predictand  classes  were  then  compared  with  each  predictor 
and  either  a  1  or  a  0  was  recorded.  The  group  means,  the  grand  means 
the  within-group  sums,  and  the  between-group  sums  were  calculated. 

The  following  values  were  obtained: 


Predictor  1:  X-q  =  0.132 

X12  =  0.097 


X.^  =  0.090 


X  =  0.106 


Predictor  2:  X ^  =  0.347 

X22=°.125 
X23  =  0.118 


X2  =  0.197 


Predictor  3:  X^-^  =  0.069 

X32  =  0.069 
X33  =  0.097 


X3  =  0.078 


SSWCXi) 

SSW(X2) 

SSW(X3) 

SSB(X1) 

SSB(X2) 

SSB(X3) 


40.9536 

63.3736 

31.2600 

0.1459 

4.8852 

0.0753 


. 
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SPWCX^)  =  13.1386 
SPW(X1X3)  =  10.5927 
SPW(X2X3)  =  13.6086 
SPB  (X1X2)  =  0.8352 

SPB(X1X3)  =  -0.0603 
SPB(X2X3)  =  -0.3312 

The  matrices  W  and  B  were  constructed  according  to  Eq .  (9)  and 

Eq .  (10)  and  the  matrix  C  =  W~^B  was  calculated  giving  the  following 
values : 

f  0.001066 
C  =  J  0.014837 
-0.008749 

Solving  the  determinantal  Eq .  (20)  for  A  gave  the  cubic 
equation : 

A3  -0.093193A2  +  0.000175A  -0.00000001  =  0  (35) 


0.005565 

0.086719 

•0.050233 


-0.000831 


-0.006205 


0.005408 


(34) 


with  the  roots: 


A  =  0 

1 

A2  =  0.091276 
A3  =  0.001917 


The  vectors  V.  representing  the  weights  on  the  predictors  are 


the  solutions  of  the  linear  homogeneous  Eq .  (21). 


To  obtain  non-trivial  solutions  to  a  set  of  these  linear  homo¬ 


geneous  equations  the  following  relation  may  be  used : 


V. 

J 


(-1) 


(36) 


• 
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where  -  determinant  of  order  2  formed  by  deleting  the  row  and 
column  containing 

The  values  of  the  weights  corresponding  to  the  three  cases  are 
listed  below: 


Case  1 


Case  2 


Case  3 


0.000157 
-0.000025 
0.000014 
x2  =  0.091276 
0.000522 
0.007738 
-0.004581 
*3  =  0.001917 
0.000040 
-0.000018 
-0.000155 


Substituting  these  values  into  Eq .  (22)  gave  the  desired 
discriminant  functions: 


Q  =  0.000157  X-  -  0.000025  X9  +  0.000014  X_  (37) 

11^3 

Q2  =  0.000522  Xl  +  0.007738  X2  -  0.004581  X3  (38) 

Q3  =  0.000040  Xx  -  0.000018  X2  -  0.000155  X3  (39) 


value 


Each  X 

i 


for  Qj  (j 


(i  =  1,2,3)  must  be  either  0  or  1 ;  therefore  the  minimum 
=  1,2,3)  is  0  and  the  maximum  value  equals: 


3 

E  V.  = 

,  l 


S 
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For  perfect  predictors,  S  would  take  on  a  value  very  close  to  1. 
Comparing  this  with  equations  (37-39)  immediately  rules  out  minimum 
temperature  and  wind  direction  as  predictors  for  smoke  emission  rate. 

For  example,  consider  Eq.  (38)  which  seems  to  be  the  most  powerful  one. 
Here  the  predictors  explain  only  about  0.47>  of  the  total  variation  in 
and  so  cannot  be  used.  The  described  procedure  can  be  said  to  be 
indirect  as  it  involves  (1)  a  calculation  of  the  smoke  emission 
rate  (Q) ,  and  (2)  an  application  of  the  calculated  Q-values  in  the 
pollution  model  to  obtain  the  level  of  smoke  concentration.  We  shall 
now  describe  a  method  that  circumvents  the  problem  of  deriving  smoke 
production  rates  and  predicts  concentrations  from  meteorological 
parameters  only. 

b)  Prediction  of  Absolute  Smoke  Concentrations  ( x  ) 

Linear  multiple  regression  was  applied  to  the  data  to  obtain 
a  prediction  equation  for  daily  smoke  concentrations.  This  method  is 
simple  and  does  not  require  vast  calculations.  Compared  with  discrimi¬ 
nant  analysis  it  has  the  weakness  of  not  accounting  for  interrelations 
between  the  variates.  Perfect  interrelationships  between  one  or  more  of 
the  variates  will  lead  to  a  singular  matrix  that  has  no  solution.  How¬ 
ever,  as  in  most  meteorological  problemsj  the  variates  are  only  partially 
interrelated,  and  so  introduce  an  error  in  the  prediction  equation  whose 
nature  is  an  over-  or  an  under-estimation,  depending  on  the  correlation 
between  the  variates. 

Variates 

Three  predictors  were  selected  to  explain  the  variation  in 
smoke  concentration  in  the  city.  Daily  values  of  these  predictors  for  six 


* 

. 

- 

* 

t 

' 


-33- 


winter  months,  October  through  March,  1967-68,  were  taken  from  Monthly 

Meteorological  Summary,  Department  of  Transport  (Canada) .  The  predict- 

and ,  valid  for  the  same  period  of  time  as  the  predictors,  was  taken 

from  Air  Pollution  Survey,  Monthly  report.  Department  of  Health.  A 

listing  of  the  variates  and  their  units  is  given  below: 

-  3  1  9 

X:  Mean  daily  concentration  of  smoke  (g.  cm.  x  10  ).  The 

initial  set  of  data  was  obtained  as  in  section  a)  of  this 
chapter . 

:  Mean  daily  wind  speed  (m.p.h.)  as  measured  10  m.  above  the 

ground.  In  the  strict  sense  the  data  are  valid  only  for  the 
Industrial  Airport  in  Edmonton. 

X 2‘.  Mean  daily  heat  island  intensity  defined  as  the  difference  in 

temperature  between  the  Industrial  and  International  Airports 
(°F) .  The  individual  values  were  obtained  from: 
l/2  (  Tmax  -  Tmin  ) 

X^:  Mean  daily  vertical  temperature  gradient  (°C(100m.)  ).  These 

standard  units  are  obtained  from  readings  at  17  m.  and  112  m. 
on  the  CN  Tower  in  downtown  Edmonton. 

Mathematical  Procedure 

We  can  write  the  regression  equation  in  the  general  form: 

Y  =  K0  +  K1X1  +  K2X2  +  K„Xn  (40) 

where  K  is  an  additive  constant  and  K-(i=l,  -  n)  are  the  weights  on 

o  1 

the  independent  variables  to  be  determined  from  the  data  in  the  sample. 
This  equation  is  the  best  linear  relation  suited  to  predict  Y  from 
Xi(i=l,  -  n)  and  is  obtained  by  the  method  of  least  squares  such  that 
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E  (Y  -  Y)  is  a  minimum.  Here  Y  are  the  observed  values.  The 
weights  can  be  computed  by  solving  simultaneously  the  normal  equations. 


NK0 

+  Kj_ 

EX1 

+  k2 

E  X2  +  ---  Kn  E  Xn  = 

E 

Y 

Ko 

E  X-i  + 

i 

K1 

2 

E  Xx 

+  K2  EXiX2  +  —  K„ 

E 

XlXn  -  2 XXY 

i  i  i 


K 

o 

IV 

■1 

+  Ki 

EXn-lXi  + 

K2  zXn-lX2  +  " 

-  ^ 

xXn-lXn 

= 

Ko 

EXn 

+ 

Kx  E 

XnXl  +  K2 

E  XnX2  + 

Kn 

x xnxn  = 

Exnx 

The  computed  sums  were  in  this  case: 


E  X1  = 

1657.4 

E  X2  =  452 

ex3  = 

77.3 

E  Xx2  = 

17530.06 

E  X22  =  2250 

-32  = 

193.55 

E  XXX2= 

3428.9 

E  XXX3=  411.91 

E  X2X3= 

282.7 

E  Y  = 

1479  E  X,Y 

=  11947.2  E  X2Y  =  4271.1 

E  X^Y  = 

1009.28 

and  the  weights: 


K  =  9.77 

o 

K  =  -0.32 

K  =  0.22 

2 

K3  =  1.68 

Substituting  the  weights  into  Eq .  (40)  gives  the  desired 
functional  relation: 

X  =  9.77  -  0.32  Xj_  +  0.22  X2  +  1.68  X3  (42) 

Data  Test 

To  test  the  stability  of  the  data  used,  six  consecutive  random 
samples  of  size  90  by  4  were  removed  from  the  original  183  by  4  data  matrix. 


. 

-35- 


A  complete  analysis  was  applied  to  each  set  with  the  following 


result : 


Table  2 


Results  of  Data  Stability  Tests 


Test 

Ko 

Ki 

k2 

k3 

r2% 

r2  % 
r  Xl/o 

r2  7 

r  X2/o 

r2  °/ 
r  x3/d 

1 

10.02 

-0.34 

0.33 

1.12 

32.4 

73.8 

8.6 

17.6 

2 

11.62 

-0.45 

0.02 

1.79 

42.3 

24.1 

0 

75.9 

3 

9.86 

-0.31 

0.06 

1.74 

38.7 

17.1 

0.3 

82.4 

4 

10.02 

-0.31 

0.28 

1.93 

50.7 

14.8 

3.6 

81.5 

5 

9.89 

-0.37 

0.30 

1.13 

41.0 

77.6 

10.2 

12.2 

6 

11.11 

-0.41 

0.24 

1.13 

31.0 

81.6 

3.9 

14.2 

Mean 

10.42 

-0.37 

0.21 

1.47 

39.4 

44.7 

4.3 

51.0 

In  column  five  is  given  the  total  variation  explained,  expressed 
in  percent.  In  column  six,  seven  and  eight  are  given  the  proportion 
of  the  variation  of  the  dependent  variable  accounted  for  by  each  independ¬ 
ent  variable,  also  in  percent. 

Judging  from  the  table  the  data  show  considerable  stability  with 
only  minor  variations  in  the  weights.  The  mean  weights  are  in  fair 
agreement  with  those  of  Eq.  (42). 

The  three  independent  variables  together  explain  about  40%  of 
the  total  variation  in  smoke  concentration.  Of  the  three,  wind  speed 
and  vertical  temperature  structure  were  of  equal  importance,  while 
the  horizontal  temperature  gradient  was  not  a  significant  predictor.  Adding 
more  variables  might  throw  some  light  on  the  somewhat  low  correlation. 

The  importance  of  relative  humidity  to  smoke  concentration  near 
the  ground  should  be  looked  into.  Improved  data  for  the  predictand 
may  also  change  the  picture.  As  of  now,  the  correlation  between  the 


, 
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the  different  sampling  stations  leaves  much  to  be  desired. 

Test 

It  was  stated  above  that  the  regression  equation  would  explain 
about  40 /o  of  the  total  variation  of  smoke.  To  test  the  stability  of 
the  equation,  an  independent  sample  comprising  one  month  of  data 
(December,  1968)  was  selected.  The  independent  variables  were  listed 
and  Eq .  (42)  solved  for  X  giving  the  daily  mean  concentration  of  smoke 

in  Edmonton.  These  values  were  then  compared  with  the  observed  values 
taken  from  Air  Pollution  Survey,  Monthly  report.  Department  of  Health 
(Figure  5  ) .  To  find  out  how  much  of  the  total  variance  the  prediction 
equation  could  explain  for  the  selected  month  ,the  differences  between 
the  predicted  and  observed  values  were  formed  with  the  mean  values  for 
the  two  sets  not  removed.  The  variance  of  the  differences  was  then 

__  o  2 

calculated.  This  came  out  to  36(yg.  m.-J)  .The  total  variation  to  be 

_  3  2 

explained  was  found  to  be  37(hg.  m.  )  .Stated  otherwise,  the  prediction 
equation  was  able  to  explain  only  about  3%  of  the  total  variation  in 
smoke  concentration.  It  was  then  suspected  that  the  difference  in  the 
means  between  the  two  sets  of  data  was  to  blame  for  the  low  value. 

In  order  to  make  the  data  independent  of  the  mean  this  was  removed 
from  each  set  of  data.  The  differences  were  again  formed  using  this 
new  set  of  data  and  the  variance  of  the  differences  was  once  more 
calculated.  This  time  the  value  came  out  to  20.0(yg.  m.  ^)2.With  the 
means  thus  removed  the  prediction  equation  explains  about  44%  of  the 
total  variation  in  smoke  concentration.  This  is  in  excellent  agreement 
with  the  value  (40%)  obtained  from  the  dependent  sample.  With  some 
reservations  because  of  the  small  test  sample  it  can  be  said  that  the 
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prediction  equation  is  quite  stable  for  daily  variations  but  fails  to 
show  skill  in  predicting  monthly  mean  smoke  concentrations. 


CHAPTER  III 


THE  HORIZONTAL  AND  VERTICAL  TEMPERATURE  DISTRIBUTION 

IN  EDMONTON  AND  CALGARY 

Introduction 

In  this  chapter  the  horizontal  and  vertical  temperature  dist¬ 
ributions  in  Edmonton  and  Calgary  are  described.  A  few  horizontal  temp¬ 
erature  traverses  were  analyzed  and  compared  with  the  results  obtained 
from  the  study  of  the  thermograph,  data. 

Many  studies  of  urban  climates  have  been  carried  out  in  the 

past  few  decades.  Outstanding  examples  are  the  fundamental  work  by 

1  2 
P.A.  Kratzer  in  his  "The  Climate  of  Cities",  and  G.A.  De  Marrais’ 

study  of  Louisville,  Kentucky.  The  latter  was  used  for  comparative 

purposes  in  the  present  study.  In  1962  W.  Klassen  did  a  comparative 

climatological  study  involving  temperature,  wind  speed  and  humidity 

3 

on  the  plain  and  in  the  river  valley  system  in  the  Edmonton  area  . 

A.  Daniels  studied  the  role  of  the  urban  heat  island  in  Edmonton 

4 

using  horizontal  temperature  traverses  . 


^P.A.  Kratzer,  The  Climate  of  Cities,  Am.  Met.  Soc.  Boston, 
Massachusets :  1956,  221  p.p. 

2 

G.A.  De  Marrais,  Vertical  Temperature  Difference  Observed 
Over  an  Urban  Area.  Am.  Met.  Society,  Bulletin,  Vol.  42,  No.  8, 
August  1961,  p.  548-554. 

W .  Klassen,  Micrometeorol ogical  Observations  in  The  North 
Saskatchewan  River  Valley  at  Edmonton,  Department  of  Transport,  Cir- 
3652  Tec-408,  May  7,  1962. 

^A.  Daniels,  The  Urban  Heat  Island  and  Air  Pollution,  Unpub¬ 
lished  M.Sc.  Thesis,  University  of  Alberta,  Edmonton,  1964,  144  p.p. 
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Thermograph  Data 

The  collection  of  thermograph  data  was  initiated  on  a  continuous 
basis  in  January  and  February,  1968.  In  the  analysis  to  follow  it  would 
have  been  preferable  to  work  with  a  full  year  of  data.  As  it  turned 
out,  however,  November  and  December,  1968,  were  never  processed  because  of 
lack  of  time  and  the  January  and  February  data  had  too  many  discontinu¬ 
ities  to  be  worth  considering.  Station  3,  located  at  St.  Stephen's 
College  on  the  campus  of  the  University  of  Alberta,  was  often  subjected 
to  vandalism  leaving  data  gaps  of  up  to  ten  to  fifteen  consecutive  days. 
Because  of  this  it  was  decided  to  omit  data  from  this  station  entirely. 
Occasional  data  gaps  were  filled  in  for  the  other  six  stations  using 
data  from  the  Industrial  Airport.  This  procedure  is  far  from  perfect 
but  will  give  a  better  estimate  of  mean  values  than  if  these  gaps 
were  left  blank.  The  author  was  now  left  with  eight  months,  March 
through  October,  1968,  of  complete  raw  data  for  six  different  thermo¬ 
graph  stations.  It  was  known  that  most  stations  were  reading  too  high 
on  sunny  days.  As  this  radiation  error  caused  by  solar  heating 

could  be  different  for  each  station  and  difficult  to  estimate  without 
extensive  calibrations, it  was  decided  to  omit  the  hours  9  a.m.  to 
5  p.m.  for  every  day.  Mechanical  errors  due  to  vandalism,  deviations 
in  the  procedure  of  changing  charts,  etc.,  were  compensated  for  by 
applying  a  correction  factor  to  the  mean  monthly  temperature  values. 

The  correction  factors  were  provided  by  Geoscience  Research  Associates, 
Limited.  These  factors,  different  for  each  month  and  each  station, 
were  obtained  by  comparing  the  raw  data  with  the  readings  from  the 
Industrial  and  International  Airports  in  periods  of  continuous  rain 
or  snow.  The  standard  deviation  of  the  differences  of  ten  hour  means 
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in  such  periods  was  ±0.4°C.  Table  3  lists  the  mean  monthly  values, 
the  corrections  and  the  corrected  mean  monthly  values  (°C.)  for  the 
six  thermograph  stations  and  the  Industrial  Airport. 

The  existence  and  intensity  of  a  heat  island  over  an  urban  area 
is  best  studied  by  measuring  the  difference  in  temperature  between  the 
city  itself  and  its  surroundings.  There  should  exist  a  well  marked 
positive  temperature  gradient  as  one  approaches  the  center  of  the  city 
from  the  outskirts.  This  difference  in  temperature  is  caused  by  dif¬ 
ferences  in  heat-storage  capacity,  differences  between  incoming  and 
outgoing  radiation  due  to  pollution,  and  the  excess  of  heat  in  the  city 
due  to  human  activities.  On  clear  calm  nights  this  difference  is  quite 
pronounced.  With  high  winds  and  cloudy  skies  it  will  decrease  accord¬ 
ingly.  To  find  out  whether  there  exists  a  significant  change  in  temper¬ 
ature  from  one  place  to  another  a  Student's  t-test  was  applied  to  the 
corrected  temperature  data. 

Statistical  Approach 

We  wanted  to  show  whether  there  exists  a  significant  change  in 
temperature  from  one  place  to  another.  If  the  Student  t-test  is  applied 
to  the  data  the  following  will  hold:  if  the  means  of  two  samples  differ 
from  each  other  so  that  t  is  greater  than  a  certain  limiting  value, 
depending  on  the  degrees  of  freedom  involved,  it  is  extremely  unlikely 
that  the  samples  are  random  samples  drawn  from  the  same  population. 

As  any  two  samples  which  are  compared  are  not  independent  random  samples 
but  are  deliberately  paired  to  reduce  all  differences  other  than  the 
ones  due  to  the  effect  being  investigated  the  normal  t-test  for  independ¬ 
ent  samples  was  ruled  out.  Instead  the  method  of  paired  samples  was 
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Obtain  the  N  differences  between  members  of  a  pair 


d .  =  X  -  Xn  . 
i  li  2i 


(1) 


where  the  subscript  1  refers  to  one  sample  and  the  subscript  2  to 

the  other.  On  the  null  hypothesis  that  the  samples  do  not  really 

differ  the  expectation  of  di  is  zero. 

2 

If  S  is  the  variance  of  the  differences  i.e. 


(2) 


the  quantity 


-  %  -1 

dN  S 


(3) 


has  the  Student-t.  distribution  with  n  degrees  of  freedom. 

Eq.  (3)  solved  for  fifteen  different  combinations  gave  the 
following  t  scores. 


5 

E . S .  Keeping ,  Introduction  to  Statistical  Inference , 
D.  van  Nostrand  Company,  Inc,,  New  York.  p.  185-186. 
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Table  4 

Results  of  Student's  t-test  for  Edmonton 


Station  Pairs  t  -  score 


2,1 

4.6 

4,1 

3.2 

5,1 

0.1 

6,1 

0.0 

7,1 

0.7 

4,2 

9.9 

5,2 

5.1 

6,2 

7.0 

7,2 

4.9 

XD,  1 

4.5 

XD,  2 

0.7 

XD,4 

15.2 

XD,  5 

8.4 

XD,  6 

9.4 

XD,7 

7.6 

Entering  tables  with  n  =  N-l  degrees  of  freedom  gave  the  limit¬ 


ing  values  of  t  of  3.5  at  the  1  °L  confidence  level  and  2.4  at  the  5 % 


level.  Thus  every  absolute  value  of  t  greater  than  these  limiting 


values  would  show  a  significant  difference  in  mean  temperature  at  these 

6 

levels  . 


Discussion 


The  t  -  score  itself  does  not  say  which  of  the  stations  com¬ 


pared  is  the  warmer.  Therefore,  the  mean  temperature  for  each  station 


together  with  its  standard  deviation  is  listed  in  Table  5. 


6 


Loc . 


cit . 
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Table.  5 

Edmonton  Thermograph  Data,  Means  and  Standard  Deviations 


Station 

Mean  °C 

SD  °  c 

1 

7.8 

5.6 

2 

8.6 

5.4 

4 

7.3 

5.4 

5 

7.8 

5.5 

6 

7.8 

5.5 

7 

7.6 

5.5 

XD  (Edmonton  Industrial 

8.7 

5.3 

Airport) 

Considering  Station  1  first  it  is  obvious  from  the  tables  that 
this  station  was  not  significantly  warmer  than  the  stations  on  the  per¬ 
iphery  in  spite  of  its  more  central  location.  The  only  exception  to 
this  is  Station  4,  located  in  the  southern  suburbs.  A  possible  ex¬ 
planation  to  the  fact  that  the  valley  station  is  not  warmer  than 
the  outer  ones  is  that  cold  air  at  nights,  and  particularly  on  calm 
and  clear  ones,  drains  down  the  valley.  Differences  in  density  would 
then  cause  a  cold  stream  of  air  to  enter  the  system  either  in  the  SW 
or  NE  portions  of  the  city.  This  stream  will  stay  cold  throughout  the 
entire  valley  as  no  air  from  the  city  is  likely  to  be  fed  in,  again  due 
to  density  differences.  It  can  be  supposed  that  this  pattern  is  dis¬ 
turbed  to  some  degree  when  overcast  conditions  are  at  hand  and  when  the 
prevailing  winds  exceed  a  certain  value.  Thus  Station  1  represents 
the  river  valley  system  which  tends  to  separate  the  city  into  at  least 
two  heat  islands.  According  to  Table  4,  Station  2,  located  at  City 
Hall  was  significantly  warmer  than  Station  1.  Unfortunately  Station  3, 
on  campus,  could  not  be  included  in  the  investigation  but  would  very 


' 
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likely  have  given  the  same  result  when  compared  with  Station  1.  It 
will  seen  from  Tables  4  and  5  that  Station  2  was  significantly 
warmer  than  any  of  the  peripheral  stations  and  this  at  the  1%  confi¬ 
dence  level.  This  was  also  true  comparing  XD  with  the  peripheral  stations. 
No  significant  difference  in  mean  temperatures  was  found  for  the 
Stations  XD  and  2  although  XD  is  situated  at  the  edge  of  a  large  field. 
Conclusion 

There  exists  a  significant  change  in  temperature  as  one  approach¬ 
es  the  center  of  the  city  froai  the  outskirts.  This  difference  is 
caused  by  differences  in  heat-storage  capacity,  differences  between 
incoming  and  outgoing  radiation  due  to  pollution,  and  the  excess  heat 
due  to  activities  in  the  city.  The  urban  heat  island  is  divided  into 
at  least  two  parts  by  the  valley  of  the  North  Saskatchewan  River. 
Temperature  Traverses 

In  order  to  supplement  the  results  from  the  analysis  of  the 
thermograph  data  a  few  temperature  traverses  were  made  over  the  city. 

A  temperature  traverse  can  be  defined  as  a  series  of  temperature  measure¬ 
ments  taken  along  a  predetermined  route.  Such  a  method  requires  an 
instrument  that  is  able  to  respond  quickly  to  small  changes  in  temperature. 
The  thermometer  used  in  the  investigation  is  described  in  Chapter  I. 

The  calculated  lag  coefficient  was  considered  small  enough  for  the 
purposes  of  the  study. 

Because  the  traverses  extended  over  several  hours  it  was  neces¬ 
sary  to  reduce  all  readings  along  the  route  to  one  constant  time,  to 
account  for  the  general  change  in  temperature  with  time.  This  was  done 


as  follows : 


m 

■\  *■ 
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To  the  measured  temperature  at  location  i  and  at  time  t  an 
amount  was  added  or  subtracted  to  give  the  temperature  at  the  constant 
time  c.  This  amount  was  calculated  from  readings  at  the  Industrial 
Airport,  and  does  not  take  into  account  the  differences  of  the  rate 
of  change  of  temperature  for  different  parts  of  the  city.  Also  the 
general  change  in  temperature  had  to  be  assumed  linear  in  time. 

Despite  these  shortcomings,  the  method  was  found  satisfactory  for  the 
purpose  of  this  study. 

Thus  stated  mathematically: 


T- 

xi,c 


•^i,t  +  ( a ,  c  "  Ta , t ) 


(4) 


where 

T.  is  the  temperature  in  °F.  at  the  constant  time  c 
1 ,  c 

Ti  t-  is  the  measured  temperature  in  °F.  at  location  i 
at  time  t 

T  is  the  Industrial  Airport  temperature  in  °F .  at 

cl  j  C 

time  c  (on  the  hour) 

Ta  t  is  the  Industrial  Airport  temperature  in  F .  at 

time  t,  determined  by  linear  interpolation  between 
two  hourly  readings.  Time  t  occurs  in  this  interval. 


In  the  following,  four  cases  are  discussed: 

Short  discussion 

Case  1.  October  10,  1968,  1115  -  1405  (Fig_. _ 6) 

Current  weather:  Wind  E  to  NW,  7  m.p.h.,  oveicast. 

this  overcast  situation  the  horizontal  temperature  distri 


* 

,  ,  -  ■  St  It  4  '■  W "!1:a 
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Figure  6 

Isotherm  Map  for  Edmonton,  October  10,  1968 
1115  -  1405  hours  MST 
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bution  showed  only  weak  gradients  in  spite  of  fairly  light  wind  speeds. 
Note  the  insignificant  influence  of  the  North  Saskatchewan  River  valley 
on  the  temperature  distribution. 

Case  2.  October  17,  1968,  0450  -  0650  (Fig.  7  ) 

Current  weather:  Wind  S,  14  m.p.h.,  clear. 

Even  this  early  morning  situation  shows  a  poorly  developed 
temperature  distribution.  The  city-country  temperature  difference 
was  not  greater  than  2°F.  The  relatively  high  wind  speed  tended  to 
make  the  distribution  of  temperature  fairly  uniform.  Note  that  the 
river  valley  once  again  had  very  little  influence  on  the  distribution. 
Case  3.  January  18,  1.969,  1815  -  2112  (Fig.  8) 

Current  weather:  Wind  NW,  2.8  m.p.h..  Variable  cloudiness. 

This  situation  shows  a  rather  well-developed  horizontal 
temperature  distribution.  Note  how  the  cold  air  drainage  in  the  river 
valley  divided  the  city  into  two  sectors.  The  temperatures  along  the 
southern  edge  were  much  lower  than  these  along  the  northern  one, 
possibly  indicating  a  lesser  city  influence  on  the  southern  suburbs. 

The  city-country  temperature  difference  was  as  large  as  8°F. 

Case  4.  January  21,  1969,  0230  -  0400  (Fig.  9  ) 

Current  weather :  Wind  SW,  2  m.p.h..  Variable  cloudiness. 

In  this  early  morning  situation  note  the  pillow  of  warm  air 
over  the  commercial  area  to  the  south  and  the  insignificant  influence 
of  the  river  valley  on  the  temperature  distribution.  Also  in  this  case 
the  city-country  temperature  difference  amounted  to  8  F. 

Vertical  Temperature  Structure 


Some  of  the  properties  of  the  vertical  temperature  structure 
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>i.  Albert 
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Figure  7 

Isotherm  Map  for  Edmonton,  October  17, 
0450  -  0650  hours  MST 


1968 


Figure  8 

Isotherm  Map  for  Edmonton,  January  18,  1969 

1815  -  2112  hours  MST 


■ 
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Figure  9 


Isotherm  Map  for  Edmonton,  January  21,  1969 
0230  -  0400  hours  MST 
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over  the  cities  of  Edmonton  and  Calgary  are  described  in  this  section. 
Two  separate  Fortran  programs  were  written  by  the  author  to  extract 
some  statistical  properties  from  the  data.  The  data  for  Edmonton 
were  collected  during  the  months  of  August,  1967,  through  July,  1968, 
and  represent  the  situation  in  the  center  of  the  city  (Fig.  1). 
Thermometers  were  mounted  at  49  feet,  186  feet  and  369  feet  on  the 
CN  Tower. 

The  Calgary  Tower  located  in  an  open  grass  covered  area  in 
the  southeastern  end  of  the  city  had  four  sensors  placed  at  10  feet, 

100  feet,  200  feet  and  300  feet  above  the  ground.  The  data  representing 
the  months  of  March,  1968,  through  September,  1968,  were  processed  simi¬ 
larly  to  the  Edmonton  data  (Fig.  4). 

Data  Processing 

The  temperature  differences  between  each  of  the  thermometers 
on  the  two  towers  were  obtained  and  set  up  as  two-dimensional  arrays 
representing  the  daily  and  hourly  values.  These  differences  were  then 
transformed  into  a  standard  unit  °C(100  m.)  ^  in  order  to  compare  them 
with  results  from  other  sites. 

The  frequency  distribution  for  each  month  using  a  temperature 
interval  of  1°C(100  m.)  ^  was  computed  and  tabulated  for  each  hour  of 
the  day. 

For  each  layer  the  monthly  mean  values  for  each  hour  were 
calculated  and  graphed.  As  previously  mentioned  in  Chapter  I,  the 
middle  thermometer  on  the  Edmonton  Tower  was  not  registering  correctly 
and  the  author  decided  to  omit  the  data  in  the  analysis.  Therefore  the 


1 
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present  discussion  is  limited  to  a  single  layer  of  thickness,  320  feet 
(  98  m.).  It  was  found  that  the  300  foot  sensor  on  the  Calgary  Tower 
was  in  error  probably  because  of  radiation  effects.  Therefore,  only  data  from 
the  lowest  three  levels  were  used  giving  a  thickness  of  190  feet  or 
approximately  58  meters  for  the  whole  structure  and  90  feet  or  approx¬ 
imately  27  meters  for  the  lowest  layer. 

Analysis 

The  analysis  is  in  two  sections:  A,  Edmonton  and  B,  Calgary. 

A.  Edmonton 

As  seen  in  Table  6,  in  which  the  average  hourly  temperature 
differences  for  alternate  hours  for  each  month  are  listed,  near- 
adiabatic.  (0.98°C/l00  m.)  conditions  prevailed  during  the  daylight 
hours  in  the  summer;  and  isothermal,  or  a  weak  positive  or  negative, 
lapse  prevailed  throughout  the  rest  of  the  year.  The  maximum  positive 
lapse  value  occurred  at  noon  ±  2  hours  due  to  maximum  incoming  solar 
radiation  at  this  time. 

Temperature  inversions  of  moderate  strength  occurred  through¬ 
out  the  year  at  night  and  in  the  early  morning  hours,  breaking  approx¬ 
imately  2-3  hours  after  sunrise  and  forming  again  at  sunset.  The  above- 
mentioned  lag  is  constant  throughout  the  year.  Note  that  because 
stabilization  began  before  sunset  there  was  no  lag  in  the  formation 
of  the  inversions  at  this  time.  The  maximum  negative  lapse  occurred 
during  the  summer  months  around  4  a.m.  and  in  winter  about  5  a.m. 

The  graphs  in  Appendix  I  reveal  an  interesting  feature,  namely  that 
during  the  summer  months,  June,  July,  August,  and  to  some  extent  Sept- 


. 


Mean  Hourly  Temperature  Gradients  in  °C(100m)  for  the  months  of  August  1967  through  July  1968  in 
Edmonton,  Alberta.  (49  ft.  to  369  ft.)  Curved  lines  through  data  indicate  sunrise  and  sunset. 
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ember  and  October,  there  is  clearly  a  secondary  maximum  around  2  p.m. 
From  a  minimum  at  noon  the  temperature  gradient  dropped  and  quickly 
approached  zero.  The  author  suggests  the  following  interpretation: 
the  forming  of  convective  clouds  around  midday  will  shield  the  earth 
enough  to  reduce  the  rate  of  heating  of  the  ground.  The  lowest  layer 
will  cool  off  somew7hat  and  a  more  homogeneous  atmosphere  will  develop 
in  the  layer  under  consideration.  To  support  this  explanation  the  total 
amount  of  cumulus,  heavy  cumulus,  cumulonimbus  and  stratocumulus  were 
plotted  for  the  hours  6  a.m.  through  5  p.m.  for  the  month  of  June, 

1968.  It  is  evident  that  these  types  of  cloud  are  at  a  maximum  between 
1  p.m.  and  3  p.m.  (Fig. 10). 

Frequency  distributions  of  the  hourly  vertical  temperature 
gradient  for  the  15-113  m.  layer  on  the  CN  Tower,  for  the  months  of 
August,  1967,  to  July,  1968,  are  shown  in  Appendix  IVa.  The  data  were 
converted  to  units  of  °C(100  m.)~^  and  grouped  in  classes  of  equal  size 
of  1°C(100  m.)-^  each.  A  pronounced  diurnal  variation,  with  frequent 
superadiabat ic  gradients  in  the  daytime,  was  evident  for  the  summer 
months  while  it  was  not  evident  for  the  winter  and  transitional  months. 
For  the  whole  year  the  most  frequent  value  of  aT  fell  in  the  class  -1 
to  0  °C(100  m.  )”■*■,  which  can  thus  be  taken  as  the  statistical  mode. 
DeMarrais  found,  in  a  study  made  in  Louisville,  Kentucky,  that,  where¬ 
as  surface  inversions  occur  regularly  over  non- urban  areas,  they  were 
comparatively  rare  over  Louisville7.  He  discovered  that  during  more 
than  half  of  the  hours  at  night,  a  lapse  rate  between  isothermal  and 

7 

DeMarrais,  op.  cit. 
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Total  amount  of  cu,  cu+,  cb  and  sc  in 
tenths  plotted  against  hours  of  the  day 
for  the  month  of  June  1968. 
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adiabatic  was  observed  in  a  150-meter  thick  layer  over  the  city.  In 
addition  he  found  that  a  lower  layer  was  considerably  less  stable  than 
a  deeper  layer. 

Except  for  the  months  August,  1967,  September,  1967,  and  July, 

1968,  it  is  obvious  that  during  about  half  of  the  hours  at  night  there 
existed,  on  the  average,  a  lapse  rate  between  isothermal  and  adiabatic. 
Taken  over  the  whole  day,  inversions  were  observed  during  approximately 
377,  of  the  total  monthly  hours  in  the  summer,  52 °/0  in  winter  and  437, 
in  the  transitional  seasons. 

Graphs  in  Appendix  II  show  the  vertical  temperature  variation  in 

the  lowest  layer.  It  was  previously  stated  that  the  middle  sensor  was 

in  error  and  thus  data  from  this  level  were  not  going  to  be  considered. 

However,  it  was  later  found  that  this  error  could  be  approximately 

corrected  by  substracting  an  average  value  of  0.7°F.  from  all  individual 

8 

readings  at  this  level  .  This  was  applied  to  data  for  the  months  of 
October-November ,  1968,  and  February-March,  1969.  The  middle  sensor  is 
located  at  a  height  of  186  feet  which  gives  the  bottom  layer  a  thickness 
of  137  feet  or  approximately  42  meters. 

The  mean  hourly  values  for  each  month  were  calculated  for  the 
corrected  sets  of  data  and  graphed  simultaneously  with  the  98-meter 
layer.  These  graphs  can  then  be  compared  with  similar  graphs  from  the 
city  of  Calgary  (Appendix  HI )  .  It  is  obvious  that  the  Calgary  daytime 
temperature  decreases  with  height,  rapidly  in  the  lower  layers  and  more 

O 

Report  of  Geoscience  Research  Associates  Limited,  "Investigation 
into  Pollution  in  Calgary  and  Edmonton",  April  1,  1967,  March  31,  1968. 
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slowly  at  greater  heights.  At  night  the  temperature  increases  with 

height  ,  the  gradient  again  being  greatest  in  the  lowest  layer.  These 

results  are  in  close  agreement  with  records  obtained  in  open  cultivated 

9 

country  .  Note  that  the  Calgary  Tower  is  located  in  the  southeastern 
part  of  the  city  in  a  valley  and  on  grass  covered  ground.  The  Edmonton 
graphs  show  a  completely  different  picture.  Here  the  curves  are  re¬ 
versed  at  night  so  that  the  lower  layer  shows  less  stability  than  the 
deeper  layer;  this  agrees  with  the  findings  of  DeMarrais  for  Louisville, 
Kentucky.  The  CN  Tower  is  located  in  central  Edmonton  and  would  more 
truly  represent  the  situation  over  the  city. 

Interpretation 

In  the  country,  and  thus  also  for  Calgary,  the  morning  sun  will 
warm  the  ground  and  the  air  immediately  above  it  much  faster  than  any 
deeper  layer  of  air.  This  will  cause  a  much  steeper  negative  temperature 
gradient  in  the  lowest  layers.  It  should  be  noted  from  the  Calgary 
graphs  that,  except  at  midday,  9T/  9t  is  very  nearly  constant  with 
respect  to  height  and  time  during  the  hours  1030  to  1400.  After  sun¬ 
set  the  lowest  layers  will  cool  off  quickly  because  of  the  greater 
loss  in  longwave  radiation.  A  stronger  inversion  will  thus  build  up 
in  the  layer  close  to  the  ground.  The  reverse  situation  at  night  in 
Edmonton  is  interpreted  by  the  author  as  a  heat  island  effect.  There 
is  apparently  enough  heat  stored  in  buildings,  streets  and  other  ob¬ 
jects  during  the  daylight  hours  to  make  up  for  a  good  part  of  the  out¬ 
going  long-wave  radiation  at  night.  Compared  with  the  country  situation 

9O.G.  Sutton,  Micrometeorology,  McGraw  Hill,  New  York,  1953, 
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this  would  of  course  develop  a  less  stable  lower  layer. 

Another  interesting  feature  of  the  Calgary  graphs  is  the  tend¬ 
ency  for  3T/  9t  to  become  independent  of  height  for  a  short  period 
around  midnight,  the  numerical  value  being  almost  the  same  as  that 
found  around  noon. 

During  any  period  of  24  hours  of  clear  weather  the  temperature 
gradient  changes  sign  twice,  shortly  after  sunrise  and  before  sunset. 

Best"-  found  that  the  evening  transition  occurred  at  much  the  same  time 
before  sunset  throughout  the  year  but  that  the  time  of  morning  transition 
varied  considerably  from  winter  to  summer,  being  later  in  the  winter  . 

The  Calgary  data  cannot  be  used  for  comparison  because  no  winter 

months  were  studied.  However,  attention  is  drawn  to  the  fact  that  the 
two  layers  reach  zero  gradient  at  the  same  time  at  both  transitions. 

The  Edmonton  graphs  (App. I)repi esenting  city  climate  seem  to  have  a 
seasonal  variation.  In  the  summer  the  morning  transition  occurred 
earlier  than  in  the  winter  but  also  the  evening  transition  had  a  season¬ 
al  variation  occurring  much  later  in  the  summer  than  in  the  winter. 

If  the  two- layer  graphs  from  Edmonton  are  compared  with  the  Calgary 
graphs  it  is  seen  that  Edmonton  shows  similarities  with  Calgary  in 
the  morning  but  that  the  evening  transition  for  the  lower  layer  in 
Edmonton  occurred  about  two  hours  later  than  that  for  the  whole  layer, 
which  is  obvious  from  the  difference  in  sunrise  -  sunset  times.  This 
is  most  likely  due  to  the  excess  of  heat  in  the  city  over  the  country. 
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B.  Calgary 

As  can  be  seen  in  Table  7,  in  which  the  average  hourly  temper¬ 
ature  differences  for  alternate  hours  for  each  month  are  listed, 
adiabatic  to  superadiabat ic  conditions  prevailed  during  the  daylight 
hours  throughout  the  seven  months.  The  lower  layer  had  the  steeper 
gradient  for  reasons  discussed  earlier.  The  maximum  positive  lapse 
value  occurred  at  about  noon  for  both  layers.  Moderate  to  strong 
inversions  occurred  throughout  the  seven  months  at  night  and  in  the 
early  morning  hours,  breaking  and  forming  around  sunrise  and  sunset 
respectively.  Again  the  gradient  in  the  lowest  layer  was  most  intense. 
The  maximum  negative  lapse  rate  occurred  in  the  summer  around  2  a.m. 
for  both  layers  and,  in  the  transitional  season,  around  3  a.m.  in  the 
lower  layer,  and  around  4  a.m.  for  the  whole  layer.  The  Calgary  graphs 
have  already  been  discussed  in  detail  in  the  previous  section. 

Frequency  distributions  of  the  hourly  vertical  temperature 
gradients,  for  the  bottom  layer  and  the  whole  layer  on  the  Calgary 
Tower,  for  the  months  March,  1968,  to  June,  1968,  were  constructed 
similar  to  those  for  Edmonton,  and  are  shown  in  Appendix  IVb .  The 
data  here  are  also  in  units  ^(lOOrn)'1  and  divided  into  classes  of  equal 
size  of  1°C  (lOOm)”1  each.  A  pronounced  diurnal  variation  was  evident 
in  both  layers.  It  was  also  evident  that  superadiabatic  conditions 
were  relatively  more  frequent  in  the  lowest  layer .  The  statistical 
mode  for  the  whole  layer  was  close  to  the  wet  adiabatic  value  of  about 
0 . 5°C (100m) “ 1  in  the  lower  layer.  No  comparison  with  the  Edmonton 
data  was  possible  because  of  differences  in  the  thicknesses  of  the 
layers. 
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Comparison  with  Data  from  Other  Sites 

Many  vertical  temperature  difference  data  obtained  from 
tower  instrumentation  in  open  country  have  been,  and  are  currently  being 

collected;  data  from  Idaho  Falls,  Idaho,  and  Rye,  Sussex,  England, 

11 

are  included  here  .  As  the  data  were  presented  in  °F  they  had  to  be 
transformed  into  standard  units  °C(100  m.)  ^  in  order  to  be  comparable. 
The  Idaho  Falls  Tower  gave  readings  only  for  a  layer  245  feet  in  depth 
or  approximately  80  meters.  The  Rye  Tower  provided  data  from  two 
layers,  the  lower  being  105  feet,  or  approximately  32  meters  thick  and 
the  entire  layer  being  300  feet  thick  or  approximately  92  meters. 
Comparison  of  Edmonton  and  Idaho  Falls,  Idaho 
The  data  in  Table  8  are  monthly  average  hourly  temperature 
differences  in  °C(100  m.)  ^  based  on  observations  from  January,  1955, 
to  May,  1958.  Superadiabatic  conditions  prevailed  throughout  most  of 
the  daylight  hours  for  the  whole  year  and  very  strong  inversions  were 
legio  at  night  and  in  the  morning  hours  throughout  the  entire  year. 

A  comparison  with  Edmonton  shows  that  the  vertical  temperature  differ¬ 
ences  over  the  city  and  over  a  typical  open  country  area  were  different. 
The  difference  in  layer  thickness  should  be  considered  but  would  in 
this  case  not  change  the  overall  picture. 

Comparison  of  Edmonton  and  Rye,  Sussex,  England 
The  data  in  Table  9  are  monthly  average  hourly  temperature 
differences  for  the  entire  layer  in  °C(100  m.)  1  based  on  observations 
from  July,  1945,  to  June,  1948.  The  positive  lapse  rates  in  the  winter 
and  transitional  seasons  were  somewhat  larger  than  those  observed  in 
Edmonton,  but  the  lapse  rates  for  the  summer  months  were  quite  similar. 

^DeMarrais,  op.  cit. 
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Source:  G.A.  DeMarrais,  "Vertical  Temperature  Difference  Observed  over  an.  Urban  Area",  Bull.  Amer 
Meteor.  Soc.,  42,  1961,  Fig.  2,  p.  549. 
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Mean  Hourly  Temperature  Gradients  in  °C(100m)_i  for  the  months  of  August  through  July  in  Rye,  Sussex, 
England.  Based  on  observations  from  July  1945  to  June  1948  (350  ft.  -  50  ft.).  Curved  lines  through 

data  indicate  sunrise  and  sunset. 
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Source:  G.A.  DeMarrais,  "Vertical  Temperature  Difference  Observed  over  an  Urban  Area",  Bull.  Amer 
Meteor.  Soc.,  42,  1961,  Fig  3,  p.  550. 
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The  intensities  of  the  inversions  at  night  and  in  the  early  morning 
were  quite  comparable.  Thus  there  was  no  pronounced  difference  be¬ 
tween  the  Rye  and  the  Edmonton  data  as  one  would  have  expected  . 

A  complicating  factor  lies  in  the  climatic  differences  between  the 
two  sites. 

Comparison  of  Calgary  and  Idaho  Falls,  Idaho 
In  spite  of  the  20  meters  difference  in  thickness  in  the  two 
layers  these  sites  have  much  in  common.  They  both  show  the  features  of 
typical  open  country  sites  with  superadiabat ic  conditions  prevailing  in 
the  daytime  and  moderate  to  strong  inversions  at  night.  Note  that 
Tables  7  and  8  show  that  the  individual  values  for  the  two  stations 
were  in  close  agreement  during  the  day  but  that  the  inversions  were 
somewhat  stronger  in  Idaho  Falls.  Note  that  Calgary  tower  is  in  a  valley. 
Comparison  of  Calgary  and  R.ye,  Sussex,  England 
In  this  case  Calgary,  more  than  Rye,  (Table  10)  showed  the 
typical  features  of  a  country  station  with  steeper  positive  lapse  rates 
in  the  daytime  and  also  steeper  negative  lapse  rates  at  night  and  in 
the  early  morning.  As  in  the  case  of  Edmonton  this  can  be  explained 
possibly  by  differences  in  climate;  Rye  to  some  extent  probably  being 
influenced  by  the  sea  which  would  at  least  reduce  the  strength  of  the 


inversions . 


< 


■ 


Mean  Hourly  Temperature  Gradients  in  °C(100m)  for  the  months  of  August  through  July  in  Rye,  Sussex, 
England.  Based  on  observations  from  July  1945  to  June  1948  (155  ft.  -  50  ft.).  Curved  lines  through 

data  indicate  sunrise  and  sunset. 
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Source:  G.A.  DeMarrais,  "Vertical  Temperature  Difference  Observed  over  an  Urban  Area",  Bull.  Amer 
Meteor.  Soc.,  42,  1961,  Fig. 7,  p.  553. 


CONCLUSIONS 


Researchers  in  the  United  States  have  obtained  useful  estimates 
of  smoke  emission  rates  for  different  cities  by  using  various  stat¬ 
istical  techniques.  However,  the  author  found  that  such  methods  were 
not  applicable  to  the  city  of  Edmonton.  M.E.  Miller  and  G.C.  Holzworth 
found  that  as  much  as  707>  of  the  total  variation  in  smoke  emission 

could  be  explained  using  mean  daily  temperature  as  a  predictor  for 
,  .  .  1 

smoke  emission  rates  .  This,  compared  with  the  author's  0.47o,  shows 
clearly  that  temperature  variations  do  not  account  for  a  significant 
portion  of  the  variation  of  smoke  output  in  Edmonton,  providing  the 
mathematical  model  from  which  the  smoke  emission  rates  were  obtained 
is  valid.  If  so,  a  reasonable  explanation  for  the  low  predictability 
is  that  space  heating,  which  is  temperature  dependent,  is  not  a  major 
source  of  smoke,  because  of  the  use  of  natural  gas.  The  slight  variation 
with  temperature  in  smoke  output  from  combustion  engines  is  not  enough 
to  be  significant.  The  failure  of  this  approach  prompted  a  test  of 
linear  multiple  regression  using  meteorological  variables  for  the 
prediction  of  mean  daily  smoke  concentrations.  The  predictors  -  mean 
daily  wind  speed,  mean  daily  heat  island  intensity  and  mean  daily 
vertical  temperature  gradient  -  were  able  to  explain  407>  of  the  total 
variation  in  mean  daily  smoke  concentrations.  Of  the  three  predictors, 
wind  speed  and  vertical  temperature  gradient  were  of  equal  importance 
while  the  horizontal  temperature  gradient  was  not  a  significant  predictor. 

A  test  with  six  random  samples  showed  the  derived  weights  on  the  independent 

^M.E.  Miller  and  G.C.  Holzworth,  An  Atmospheric  Diffusion  Model 
for  Metropolitan  Areas.  Air  Resources  Field  Research  Office,  ESSA 
Cincinnati,  Ohio,  1967,  p.p.  46  -  50. 
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variables  to  be  very  stable.  A  second  test  on  an  independent  data 
sample  showed  that  the  prediction  equation  was  quite  stable  for  daily 
variations  but  that  it  failed  to  show  skill  in  predicting  monthly  mean 
smoke  concentrations. 

A  t-test,  using  the  method  of  paired  samples  applied  to  the 
Edmonton  thermograph  data,  showed  the  existence  of  a  significant 
change  in  temperature  as  one  approached  the  center  of  the  city  from 
the  outskirts.  This  difference  was  attributed  to  differences  in  heat- 
storage  capacity,  differences  between  incoming  and  outgoing  radiation 
because  of  pollution,  and  the  excess  heat  due  to  activities  in  the 
city.  The  urban  heat  island  was  divided  into  at  least  two  parts  by 
the  river  valley.  The  results  obtained  from  the  thermograph  data  were 
supplemented  in  a  few  instances  by  horizontal  temperature  traverses. 

Analysis  of  the  vertical  temperature  distribution  in  Edmonton 
revealed  that  near-adiabatic  conditions  prevailed  during  the  daylight 
hours  in  the  summer;  and  isothermal,  or  weak  positive  or  negative 
lapse  prevailed  throughout  the  rest  of  the  year.  Temperature  inversions 
of  moderate  strengths  occurred  throughout  the  year  at  night  and  in  the 
early  morning  hours. 

The  same  analysis  applied  to  Calgary  showed  that  adiabatic  to 
superadiabat ic  conditions  prevailed  during  the  daylight  hours  through¬ 
out  the  seven  months  that  were  studied.  Moderate  to  strong  inversions 
occurred  throughout  these  months  at  night  and  in  the  early  morning  hours. 
The  results  of  the  vertical  temperature  distribution  from  Edmonton  and 
Calgary  were  compared  with  similar  results  from  other  sites  and  showed 
that  Edmonton,  more  so  than  Calgary,  represented  the  typical  city 
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situation. 

peripheral 


A  possible  explanation  for  this  difference  may  be  the 
location  of  the  Calgary  instrument  tower. 
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APPENDIX  I 


Graphs  of  the  hourly  vertical  temperature  gradient  in 
°C(100  m.)  ^  between  15  m.  and  113  m.  on  the  CN  Tower  in  central 
Edmonton.  Horizontal  dashed  line  marks  the  dry  adiabatic  lapse 
rate.  Vertical  bars  on  the  abscissa  mark  the  times  of  sunrise 


and  sunset. 
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APPENDIX  II 


Graphs  of  the  hourly  vertical  temperature  gradient  in 
°C(100  m.)  ^  for  the  layers  15  m,  to  113  m.  (crosses)  and  15  m. 
to  42  m.  (circles)  on  the  CN  Tower  in  central  Edmonton.  Hor¬ 
izontal  dashed  line  marks  the  dry  adiabatic  lapse  rate.  Vertical 
bars  on  the  abscissa  mark  the  times  of  sunrise  and  sunset. 
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APPENDIX  III 


Graphs  of  the  hourly  vertical  temperature  gradient  in 
°C(100  m.)  ^  for  the  layers  3  m.  to  61  m.  (crosses)  and  3  m. 
31  m.  (circles)  on  the  Calgary  Tower  in  Southeastern  Calgary. 
Horizontal  dashed  line  marks  the  dry  adiabatic  lapse  rate. 
Vertical  bars  on  the  abscissa  mark  the  times  of  sunrise  and 


sunset . 
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APPENDIX  IV  a 


Frequency  distributions  of  vertical  temperature 
gradients  between  15  m.  and  113  m.  on  the  CN  Tower  in  central 
Edmonton  for  the  months  August  1967  through  July  1968. 

All  temperature  differences  are  converted  to  a  standard 
height  interval  of  100  m. 

Note  that  the  typed  numbers  are  frequencies  for  1°C 
temperature-gradient  intervals.  The  frequencies  listed  for 
-0 . 99 °C (100m)  \  which  indicate  the  dry  adiabatic  lapse  rate, 
have  to  be  added  to  the  row  below  to  give  the  frequencies  in  the 
temperature  interval  -0.99  to  0.01. 
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APPENDIX  IV  b 


Frequency  distributions  of  vertical  temperature  gradients 
between  3  m.  and  61  m.  and  3  m.  and  31  m.  on  the  Calgary  Tower  in 
southeastern  Calgary.  All  temperature  differences  are  converted  to 
a  standard  height  interval  of  100  m. 

Note  that  the  typed  numbers  are  frequencies  for  1°C 
temperature-gradient  intervals.  The  frequencies  listed  for  -0.99°C (100m) 
which  indicate  the  dry  adiabatic  lapse  rate,  have  to  be  added  to  the  row 
below  to  give  the  frequencies  in  the  temperature  interval  -0.99  to  0.01. 
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